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Observational cosmology is in its "golden age" with a vast amount of recent data on the distribu- 
tion of matter and light in the universe. This data can be used to probe theories of the very early 
universe. It is small amplitude cosmological fluctuations which encode the information about the 
very early universe and relate it to current data. Hence, a central topic in these lectures is the "the- 
ory of cosmological perturbations", the theory which describes the generation of inhomogeneities 
in the very early universe and their evolution until the current time. I will apply this theory to 
three classes of models of the very early universe. The first is "Inflationary Cosmology", the cur- 
rent paradigm for understanding the early evolution of the universe. I will review the successes 
of inflationary cosmology, but will also focus on some conceptual challenges which inflationary 
cosmology is facing, challenges which motivate the search for possible alternatives. I will intro- 
duce two alternative scenarios, the "Matter Bounce" model and "String Gas Cosmology", and I 
will discuss how cosmological fluctuations which can explain the current data are generated in 
those models. 
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1. Introduction 

1.1 Overview 

Observational cosmology is currently in its "Golden Age". A wealth of new observational 
results are being uncovered. The cosmic microwave background (CMB) has been measured to 
high precision, the distribution of visible and of dark matter is being mapped out to greater and 
greater depths. New observational windows to probe the structure of the universe are opening up. 
To explain these observational results it is necessary to consider processes which happened in the 
very early universe. 

The inflationary scenario [Q] (see also [|[ |], ||]) is the current paradigm for the evolution of the 
very early universe. Inflation can explain some of the puzzles which the previous paradigm - Stan- 
dard Big Bang cosmology - could not address. More importantly, however, inflationary cosmology 
gave rise to the first explanation for the origin of inhomogeneities in the universe based on causal 
physics [§] (see also [0, ||, §])). The theory of structure formation in inflationary cosmology pre- 
dicted the detailed shape of the angular power spectrum of CMB anisotropics, a prediction which 
was verified many years later observationally 

In spite of this phenomenological success, inflationary cosmology is not without its concep- 
tual problems. These problems motivate the search for alternative proposals for the evolution of the 
early universe and for the generation of structure. These alternatives must be consistent with the 
current observations, and they must make predictions with which they can be observationally dis- 
tinguished from inflationary cosmology. In these lectures I will discuss two alternative scenarios, 
the "Matter Bounce" (see e.g. [[Hj]) and "String Gas Cosmology" [11]. 



I begin these lectures with a brief survey of the Standard Big Bang (SBB) model and its 
problems. Possibly the most important problem is the absence of a structure formation scenario 
based on causal physics. I will then introduce inflationary cosmology and the two alternatives and 
show how they address the conceptual problems of the SBB model, in particular how structure is 
formed in these models. 

Since the information about the early universe is encoded in the spectrum of cosmological fluc- 
tuations about the expanding background, it is these fluctuations which provide a link between the 
physics in the very early universe and current cosmological observations. Section 2 is devoted to an 
overview of the theory of linearized cosmological perturbations, the main technical tool in modern 
cosmology. The analysis of this section is applicable to all early universe theories. In Section 3 
I give an overview of inflationary cosmology, focusing on the basic principles and emphasizing 
recent progress and conceptual problems. 

The conceptual problems of inflationary cosmology motivate the search for alternatives. In 
Section 4 the "Matter Bounce" alternative is discussed, in Section 5 the "String Gas" alternative. 

1.2 Standard Big Bang Model 

Standard Big Bang cosmology is based on three principles. The first is the assumption that 
space is homogeneous and isotropic on large scales. The second is the assumption that the dynamics 
of space-time is described by the Einstein field equations. The third basis of the theory is that matter 
can be described as a superposition of two classical perfect fluids, a radiation fluid with relativistic 
equation of state and pressure-less (cold) matter. 
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The first principle of Standard Cosmology implies that the metric takes the following most 
general form for a homogeneous and isotropic four-dimensional universe 

ds 2 = dt 2 -a{tfdx 2 , (1.1) 

where t is physical time, x denote the three co-moving spatial coordinates (points at rest in an 
expanding space have constant co-moving coordinates), and the scale factor a(t) is proportional to 
the radiius of space. For simplicity we have assumed that the universe is spatially flat. 
The expansion rate H (t) of the universe is given by 

H(t) = -, (1.2) 
a 

where the overdot represents the derivative with respect to time. The cosmological redshift z(t) 
at time t yields the amount of expansion which space has undergone between the time t and the 
present time t : 

z(t) + l = ^. (1-3) 
a(t) 

In the case of a homogeneous and isotropic metric ( |1.1[ ), the Einstein equations which describe 
how matter induces curvature of space-time reduce to a set of ordinary differential equations, the 
Friedmann equations. Written for simplicity in the case of no cosmological constant and no spatial 
curvature they are 

H 2 = —p, (1.4) 

where p is the energy density of matter and G is Newton's gravitational constant, and 

a _ AnG 
a 

or equivalently 



{p + 3p) (1.5) 



p = -3H(p+p), (1.6) 

where p is the pressure density of matter. 

The third principle of Standard Cosmology says that the energy density and pressure are the 
sums of the contributions from cold matter (symbols with subscript "m") and radiation (subscripts 
V): 

P = Pm + Pr (1-7) 

P = Pm+Pr 

with p m = and p r = l/3p r . 

The first principle of Standard Cosmology, the homogeneity and isotropy of space on large 
scales, was initially introduced soleley to simplify the mathematics. However, it is now very well 
confirmed observationally by the near isotropy of the cosmic microwave background and by the 
convergence to homogeneity in the distribution of galaxies as the length scale on which the universe 
is probed increases. 

Standard Cosmology rests on three observational pillars. Firstly, it explains Hubble's redshift- 
distance relationship. Secondly, and most importantly for the theme of these lectures, it predicted 
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Figure 1: A sketch of the temperature (vertical axis) - time (horizontal axis) relation in Standard Cosmology. 
As the beginning of time is approached, the temperature diverges. 

the existence and black body nature of the cosmic microwave background (CMB). The argument is 
as follows. If we consider regular matter and go back in time, the temperature of matter increases. 
Eventually, it exceeds the ionization temperature of atoms. Before that time, matter was a plasma, 
and space was permeated by a thermal gas of photons. A gas of photons will thus remain at the 
current time. The photons last scatter at a time t rec , the "time of recombination" which occurs at 
a redshift of about 10 3 . After that, the gas of photons remains in a thermal distribution with a 
temperature T which redshifts as the universe expands, i.e. T ~ a -1 . The CMB is this remnant 
gas of photons from the early universe. The precision measurement of the black body nature of 



the CMB [12] can be viewed as the beginning of the "Golden Age" of observational cosmology. 
The third observational pillar of Standard Cosmology is the good agreement between the predicted 
abundances of light elements and the observed ones. This agreement tells us that Standard Cos- 
mology is a very good description of how the universe evolved back to the time of nucleosynthesis. 
Any modifications to the time evolution of Standard Cosmology must take place before then. 

At the present time to matter is dominated by the cold pressureless component. Thus, as we 
go back in time, the universe is (except close to the present time when it appears that a residual 
cosmological constant is beginning to dominate the cosmological dynamics) initially in a matter- 
dominated phase during which a(t) ~ ? 2 / 3 . Since the energy density in cold matter scales as a~ 3 
whereas that of relativistic radiation scales as a~ 4 ', there will be a time before which the energy 
density in radiation was larger than that of cold matter. The time when the two energy densities are 
equal is t eq and corresponds to a redshift of around z = 10 . For t < t eq the universe is radiation- 
dominated and a(t) ~ t 1 / 2 . As we go back into the past the density and temperature increase 
without bound and a singularity is reached at which energy density, temperature and curvature are 
all infinite. This is the "Big Bang" (see Fig. [l]for a sketch of the temperature/time relationship in 
Standard Cosmology). 

1.3 Problems of the Standard Big Bang 

Obviously, the assumptions on which Standard Cosmology is based break down long before 
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Figure 2: Sketch illustrating the horizon problem of Standard Cosmology: our past horizon at t rec is larger 
than the causal horizon (forward light cone) at that time, thus making it impossible to causally explain the 
observed isotropy of the CMB. 

the singularity is reached. This is the "singularity problem" of the model, a problem which also 
arises in inflationary cosmology. It is wrong to say that Standard Cosmology predicts a Big Bang. 
Instead, one should say that Standard Cosmology is incomplete and one does not understand the 
earliest moments of the evolution of the universe. 

There are, however, more "practical" "problems" of the Standard Big Bang model - problems 
in the sense that the model is unable to explain certain key features of the observed universe. The 
first such problem is the "horizon problem" within Standard Cosmology there is no possi- 
ble explanation for the observed homogeneity and isotropy of the CMB. Let us consider photons 
reaching us from opposite angles in the sky. As sketched in Figure ||, the source points for these 
photons on the last scattering surface are separated by a distance greater than the horizon at that 
time. Hence, there is no possible causal mechanism which can relate the temperatures at the two 
points. 

A related problem is the "size problem": if the spatial sections of the universe are finite, 
then the only length scale available at the Planck time t p i is the Planck length. However, it we 
extrapolate the size of our currently observed horizon back to when the temperature was equal 
to the Planck temperature, then the corresponding wavelength is larger than 1/xm, many orders 
of magnitude larger than the Planck length. Standard Cosmology offers no explanation for these 
initial conditions. 

The third mystery of Standard Cosmology concerns the observed degree of spatial flatness of 
the universe. At the current time the observed energy density equals the "critical" energy density 
p c of a spatially flat universe to within 10% or better. However, in Standard Cosmology p = p c is 
an unstable critical point in an expanding universe. As the universe expands, the relative difference 
between p and p c increases. This can be seen by taking the Friedmann equation in the presence of 
spatial curvature 



H 2 + eT 2 



(1.8) 
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Figure 3: Sketch illustrating the formation of structure problem of Standard Cosmology: the physical wave- 
length of fluctuations on fixed comoving scales which correspond to the large-scale structures observed in 
the universe is larger than the horizon at the time t eq of equal matter and radiation, the time when matter 
fluctuations begin to grow. Hence, it is impossible to explain the origin of non-trivial correlations of the 
seeds for the fluctuations which had to have been present at that time. 



where 



{aTf 



(1.9) 



T being the temperature and k the curvature constant which is k = ± 1 or k = for closed, open 



or flat spatial sections, and comparing (1.8) with the corresponding equation in a spatially flat 
universe (k = and p c replacing p . If entropy is conserved (as it is in Standard Cosmology) then e 
is constant and we obtain 



P-Pc 



3 ST 2 



:~ 2 . (i.io) 

p c SkG p c 

Hence, to explain the currently observed degree of spatial flatness, the initial spatial curvature had 
to have been tuned to a very high accuracy. This is the "flatness problem". 

We observe highly non-random correlations in the distribution of galaxies in the universe. The 
only force which can act on the relevant distances is gravity. Gravity is a weak force, and therefore 
the seed fluctuations which develop into the observed structures had to have been non-randomly 
distributed at the time t eq , the time when gravitational clustering begins (see Section 2). However, 
as illustrated in Figure ||, the physical length corresponding to the fluctuations which we observe 
today on the largest scales (they have constant comoving scale) is larger than the horizon at that 
time. Hence, there can be no causal generation mechanism for these perturbations 1 . This is the 
"fluctuation problem" of Standard Cosmology. 

There are also more conceptual problems: Standard Cosmology is based on treating matter 
as a set of perfect fluids. However, we know that at high energies and temperatures a classical 
description of matter breaks down. Thus, Standard Cosmology must break down at sufficiently 
high energies. It does not contain the adequate matter physics to describe the very early uni- 
verse. Similarly, the singularity of space-time which Standard Cosmology contains corresponds 



This is the usual textbook argument. The student is invited to find (at least) two flaws in this argument. 
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Figure 4: A sketch showing the time line of inflationary cosmology. The period of accelerated expan- 
sion begins at time U and end at tR. The time evolution after tR corresponds to what happens in Standard 
Cosmology. 



to a breakdown of the assumptions on which General Relativity is based. This is the "singularity 
problem" of Standard Cosmology. 

All of the early universe scenarios which I will discuss in the following provide solutions 
to the formation of structure problem. They can successfully explain the wealth of observational 
data on the distribution of matter in the universe and on the CMB anisotropics. The rest of these 
lectures will focus on this point. However, the reader should also ask under which circumstances 
the scenarios to be discussed below address the other problems mentioned above. 

1.4 Inflation as a Solution 

The idea of inflationary cosmology is to add a period to the evolution of the very early universe 
during which the scale factor undergoes accelerated expansion - most often nearly exponential 
growth. The time line of inflationary cosmology is sketched in Figure |j. The time ?, is the beginning 
of the inflationary period, and tR is its end (the meaning of the subscript R will become clear later). 
Although inflation is usually associated with physics at very high energy scales, e.g. E ~ 10 16 Gev, 
all that is required from the initial basic considerations is that inflation ends before the time of 
nucleosynthesis. 

During the period of inflation, the density of any pre-existing particles is red-shifted. Hence, if 
inflation is to be viable, it must contain a mechanism to heat the universe at tR, a "reheating" mech- 
anism - hence the subscript R on tR. This mechanism must involve dramatic entropy generation. It 
is this non-adiabatic evolution which leads to a solution of the flatness problem, as the reader can 
verify by inspecting equation ( 1.10| ) and allowing for entropy generation at the time tR. 



A space-time sketch of inflationary cosmology is given in Figure ^| . The vertical axis is time, 
the horizontal axis corresponds to physical distance. Three different distance scales are shown. The 
solid line labelled by k is the physical length corresponding to a fixed comoving perturbation. The 
second solid line (blue) is the Hubble radius 



l H {t)=H-\t). (1.11) 

The Hubble radius separates scales where microphysics dominates over gravity (sub-Hubble scales) 
from ones on which the effects of microphysics are negligible (super-Hubble scales) 2 . Hence, a 
necessary requirement for a causal theory of structure formation is that scales we observe today 
originate at sub-Hubble lengths in the early universe. The third length is the "horizon", the forward 



2 This statement will be demonstrated later in these lectures. 
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Figure 5: Space-time sketch of inflationary cosmology. The vertical axis is time, the horizontal axis corre- 
sponds to physical distance. The solid line labelled k is the physical length of a fixed comoving fluctuation 
scale. The role of the Hubble radius and the horizon are discussed in the text. 

light cone of our position at the Big Bang. The horizon is the causality limit. Note that because 
of the exponential expansion of space during inflation, the horizon is exponentially larger than the 
Hubble radius. It is important not to confuse these two scales. Hubble radius and horizon are the 
same in Standard Cosmology, but in all three early universe scenarios which will be discussed in 
the following they are different. 

From Fig. |5] it is clear that provided that the period of inflation is sufficiently long, all scales 
which are currently observed originate as sub-Hubble scales at the beginning of the inflationary 
phase. Thus, in inflationary cosmology it is possible to have a causal generation mechanism of 
fluctuations [7|, g]. Since matter pre-existing at ?, is redshifted away, we are left with a matter 
vacuum. The inflationary universe scenario of structure formation is based on the hypothesis that 
all current structure originated as quantum vacuum fluctuations. From Figure Q it is also clear that 
the horizon problem of standard cosmology can be solved provided that the period of inflation lasts 
sufficiently long. The reader should convince him/herself that the required period of inflation is 
about 50H~ l if inflation takes place at an enegy scale of about 10 16 GeV. Inflation thus solve both 
the horizon and the structure formation problems. 

To obtain exponential expansion of space in the context of Einstein gravity the energy density 
must be constant. Thus, during inflation the total energy and size of the universe both increase ex- 
ponentially. In this way, inflation can solve the size and entropy problems of Standard Cosmology. 

To summarize the main point concerning the generation of cosmological fluctuations (the main 
theme of these lectures) in inflationary cosmology: the first crucial criterium which must be satis- 
fied in order to have a successful theory of structure formation is that fluctuation scales originate 
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inside the Hubble radius. In inflationary cosmology it is the accelerated expansion of space during 
the inflationary phase which provides this success. In the following we will emphasize what is 
responsible for the corresponding success in the two alternative scenarios which we will discuss. 

1.5 Matter Bounce as a Solution 

The first alternative to cosmological inflation as a theory of structure formation is the "matter 
bounce" , an alternative which is not yet well appreciated (for an overview the reader is referred to 



[10]). The scenario is based on a cosmological background in which the scale factor a(t) bounces 
in a non-singular manner. 

Figure ^ shows a space-time sketch of such a bouncing cosmology. Without loss of generality 
we can adjust the time axis such that the bounce point (minimal value of the scale factor) occurs at 
t = 0. There are three phases in such a non-singular bounce: the initial contracting phase during 
which the Hubble radius is decreasing linearly in \t\, a bounce phase during which a transition 
from contraction to expansion takes place, and thirdly the usual expanding phase of Standard Cos- 
mology. There is no prolonged inflationary phase after the bounce, nor is there a time-symmetric 
deflationary contracting period before the bounce point. As is obvious from the Figure, scales 
which we observe today started out early in the contracting phase at sub-Hubble lengths. The 
matter bounce scenario assumes that the contracting phase is matter-dominated at the times when 
scales we observe today exit the Hubble radius. A model in which the contracting phase is the 
time reverse of our current expanding phase would obey this condition. The assumption of an 
initial matter-dominated phase will be seen later in these lectures to be important if we want to 
obtain a scale-invariant spectrum of cosmological perturbations from initial vacuum fluctuations 



|[13U!4 US- 
Let us make a first comparison with inflation. A non-deflationary contracting phase replaces 
the accelerated expanding phase as a mechanism to bring fixed comoving scales within the Hubble 
radius as we go back in time, allowing us to consider the possibility of a causal generation mech- 
anism of fluctuations. Starting with vacuum fluctuations, a matter-dominated contracting phase 
is required in order to obtain a scale-invariant spectrum. This corresponds to the requirement in 
inflationary cosmology that the accelerated expansion be nearly exponential. 

With Einstein gravity and matter satisfying the usual energy conditions it is not possible to 
obtain a non-singular bounce. However, as mentioned before, it is unreasonable to expect that 
Einstein gravity will provide a good description of the physics near the bounce. There are a large 
number of ways to obtain a non-singular bouncing cosmology. To mention but a few, it has been 
shown that a bouncing cosmology results naturally from the special ghost-free higher derivative 



gravity Lagrangian introduced in Q16||. Bounces also arise in the higher-derivative non-singular 



universe construction of in "mirage cosmology" (see e.g. [18]) which is the cosmology on 
a brane moving through a curved higher-dimensional bulk space-time (the time dependence for a 
brane observer is induced by the motion through the bulk), or - within the context of Einstein gravity 
- by making use of "quintom matter" (matter consisting of two components, one with regular kinetic 
term in the action, the other one with opposite sign kinetic action) JT9[]. For an in-depth review of 
ways of obtaining bouncing cosmologies see [2C]. Very recently, it has been shown [ p] ] that 



a bouncing cosmology can easily emerge from Horava-Lifshitz gravity [22], a new approach to 
quantizing gravity. 
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Figure 6: Space-time sketch in the matter bounce scenario. The vertical axis is conformal time 77, the 
horizontal axis denotes a co-moving space coordinate. Also, _1 denotes the co-moving Hubble radius. 



In the matter bounce scenario the universe begins cold and therefore large. Thus, the size 
problem of Standard Cosmology does not arise. As is obvious from Figure ^ there is no horizon 
problem for the matter bounce scenario as long as the contracting period is long (to be specific, of 
similar duration as the post-bounce expanding phase until the present time). By the same argument, 
it is possible to have a causal mechanism for generating the primordial cosmological perturbations 
which evolve into the structures we observe today. Specifically, as will be discussed in the section 
of the matter bounce scenario, if the fluctuations originate as vacuum perturbations on sub-Hubble 
scales in the contracting phase, then the resulting spectrum at late times for scales exiting the 



Hubble radius in the matter-dominated phase of contraction is scale-invariant [13, 14, 15]. The 
propagation of infrared (IR) fluctuations through the non-singular bounce was analyzed in the case 
of the higher derivative gravity model of [16] in [23], in mirage cosmology in [pf], in the case 



of the quintom bounce in [24, 25] and for a Horava-Lifshitz bounce in [26]. The result of these 



studies is that the scale-invariance of the spectrum before the bounce goes persists after the bounce 
as long as the time period of the bounce phase is short compared to the wavelengths of the modes 
being considered. Note that if the fluctuations have a thermal origin, then the condition on the 



background cosmology to yield scale-invariance of the spectrum of fluctuations is different [27] 



1.6 String Gas Cosmology as a Solution 



String gas cosmology [f2f|] (see also p9(], and see pTT| ] for a comprehensive review) is a toy 
model of cosmology which results from coupling a gas of fundamental strings to a background 
space-time metric. It is assumed that the spatial sections are compact. It is argued that the universe 
starts in a quasi-static phase during which the temperature of the string gas hovers at the Hagedorn 
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Figure 7: The dynamics of string gas cosmology. The vertical axis represents the scale factor of the universe, 
the horizontal axis is time. Along the horizontal axis, the approximate equation of state is also indicated. 
During the Hagedorn phase the pressure is negligible due to the cancellation between the positive pressure 
of the momentum modes and the negative pressure of the winding modes, after time tR the equation of state 
is that of a radiation-dominated universe. 



value j30|], the maximal temperature of a gas of closed strings in thermal equilibrium. The string gas 
in this early phase is dominated by strings winding the compact spatial sections. The annihilation 
of winding strings will produce string loops and lead to a transition from the early quasi-static 
phase to the radiation phase of Standard Cosmology. Fig. [7] shows a sketch of the evolution of the 
scale factor in string gas cosmology. 

In Figure || we sketch the space-time diagram in string gas cosmology. Since the early Hage- 
dorn phase is quasi-static, the Hubble radius is infinite. For the same reason, the physical wave- 
length of fluctuations remains constant in this phase. At the end of the Hagedorn phase, the Hubble 
radius decreases to a microscopic value and makes a transition to its evolution in Standard Cosmol- 
ogy- 

Once again, we see that fluctuations originate on sub-Hubble scales. In string gas cosmology, 
it is the existence of a quasi-static phase which leads to this result. As will be discussed in the 
section on string gas cosmology, the source of perturbations in string gas cosmology is thermal: 
string thermodynamical fluctuations in a compact space with stable winding modes in fact leads to 
a scale-invariant spectrum [[31}] . 



2. Theory of Cosmological Perturbations 

The key tool which is used in modern cosmology to connect theories of the very early universe 
with cosmological observations today is the theory of cosmological perturbations. In the following, 



we will give an overview of this theory (similar to the overview in [33] which in turn is based on the 
comprehensive review [0]). We begin with the analysis of perturbations in Newtonian cosmology, 
a useful exercise to develop intuition and notations. 
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Figure 8: Space-time diagram (sketch) showing the evolution of fixed co-moving scales in string gas cos- 
mology. The vertical axis is time, the horizontal axis is physical distance. The solid curve represents the 
Einstein frame Hubble radius H 1 which shrinks abruptly to a micro-physical scale at tR and then increases 
linearly in time for t > tR. Fixed co-moving scales (the dotted lines labeled by k\ and ko) which are currently 
probed in cosmological observations have wavelengths which are smaller than the Hubble radius before tR. 
They exit the Hubble radius at times tj(k) just prior to tR, and propagate with a wavelength larger than the 
Hubble radius until they reenter the Hubble radius at times tf(k). 



2.1 Newtonian Theory 

The growth of density fluctuations is a consequence of the purely attractive nature of the grav- 
itational force. Imagine (first in a non-expanding background) a density excess 8p localized about 
some point x in space. This fluctuation produces an attractive force which pulls the surrounding 
matter towards x. The magnitude of this force is proportional to 8p. Hence, by Newton's second 
law 

8p~G8p, (2.1) 

where G is Newton's gravitational constant. Hence, there is an exponential instability of flat space- 
time to the development of fluctuations. 

Obviously, in General Relativity it is inconsistent to consider density fluctuations in a non- 
expanding background. If we consider density fluctuations in an expanding background, then the 



expansion of space leads to a friction term in (2.1). Hence, instead of an exponential instability 



to the development of fluctuations, fluctuations in an expanding Universe will grow as a power 
of time. We will now determine what this power is and how it depends both on the background 
cosmological expansion rate and on the length scale of the fluctuations. 
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We first consider the evolution of hydrodynamical matter fluctuations in a fixed non-expanding 
background. In this context, matter is described by a perfect fluid, and gravity by the Newtonian 
gravitational potential (p. The fluid variables are the energy density p, the pressure p, the fluid 
velocity v, and the entropy density s. The basic hydrodynamical equations are 

p + V-(pv) = 
1 

v + (v • V)v + —Vp + V<p = 

V 2 <p = \%Gp (2.2) 
s+(\-V)s = 

P = P{p,s). 

The first equation is the continuity equation, the second is the Euler (force) equation, the third 
is the Poisson equation of Newtonian gravity, the fourth expresses entropy conservation, and the 
last describes the equation of state of matter. The derivative with respect to time is denoted by an 
over-dot. 

The background is given by the energy density p Q , the pressure p , vanishing velocity, constant 
gravitational potential <po and constant entropy density sq. As mentioned above, it does not satisfy 
the background Poisson equation. 

The equations for cosmological perturbations are obtained by perturbing the fluid variables 
about the background, 

p = po + Sp 
v = <5v 

p = po + Sp (2.3) 
<p = (p + S(p 

s = s + 8s, 

where the fluctuating fields 8p ,S\,8p,S(p and 8s are functions of space and time, by inserting 



these expressions into the basic hydrodynamical equations (2.2), by linearizing, and by combining 
the resulting equations which are of first order in time. We get the following differential equations 
for the energy density fluctuation 8p and the entropy perturbation 8s 

Sp - c 2 V 2 8p - 47iGp 8p = aV 2 8s (2.4) 
8s = 0, 

where the variables c 2 and a describe the equation of state 

Sp = c 2 8p + 08S (2.5) 

with 

r 2 = I 

Sp 



d = (S?) L (2-6) 



denoting the square of the speed of sound. 
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The fluctuations can be classified as follows: If 8s vanishes, we have adiabatic fluctuations. 
If 8s is non- vanishing but 8p = 0, we speak of an entropy fluctuation. 

The first conclusions we can draw from the basic perturbation equations (2A) are that 

1) entropy fluctuations do not grow, 

2) adiabatic fluctuations are time-dependent, and 

3) entropy fluctuations seed adiabatic ones. 

All of these conclusions will remain valid in the relativistic theory. 

Since the equations are linear, we can work in Fourier space. Each Fourier component 8pk(t) 
of the fluctuation field 8p(x,t) evolves independently. In the case of adiabatic fluctuations, the 
cosmological perturbations are described by a single field which obeys a second order differential 
equation and hence has two fundamental solutions. We will see that this conclusion remains true 
in the relativistic theory. 

Taking a closer look at the equation of motion for 8p, we see that the third term on the left 
hand side represents the force due to gravity, a purely attractive force yielding an instability of flat 
space-time to the development of density fluctuations (as discussed earlier, see (2T)). The second 
term on the left hand side of (^4j) represents a force due to the fluid pressure which tends to set 
up pressure waves. In the absence of entropy fluctuations, the evolution of Sp is governed by the 
combined action of both pressure and gravitational forces. 

Restricting our attention to adiabatic fluctuations, we see from (2.4) that there is a critical 
wavelength, the Jeans length, whose wavenumber kj is given by 

kj = ( 4*Gpo )1/2 (2?) 

Fluctuations with wavelength longer than the Jeans length (k <^kj) grow exponentially 

8p k {t) ~ e Wkt with (O k ~ 4(7iGp ) 1/2 (2.8) 

whereas short wavelength modes (k S> kj) oscillate with frequency (Ou ~ c s k. Note that the value of 
the Jeans length depends on the equation of state of the background. For a background dominated 
by relativistic radiation, the Jeans length is large (of the order of the Hubble radius H~ l (t)), whereas 
for pressure-less matter it vanishes. 

Let us now improve on the previous analysis and study Newtonian cosmological fluctuations 
about an expanding background. In this case, the background equations are consistent (the non- 
vanishing average energy density leads to cosmological expansion). However, we are still neglect- 
ing general relativistic effects (the fluctuations of the metric). Such effects turn out to be dominant 
on length scales larger than the Hubble radius H~ l (t), and thus the analysis of this section is appli- 
cable only to smaller scales. 

The background cosmological model is given by the energy density po(t), the pressure po(t), 
and the recessional velocity vo = H(t)x, where x is the Euclidean spatial coordinate vector ("phys- 
ical coordinates"). The space- and time-dependent fluctuating fields are defined as in the previous 
section: 

p(t,x) = p (t)(l + 8 e {t,x)) 

\(t,x) = \ (t,x) + 8\{t,x) (2.9) 
p(t,x) = p (t) + 8p{t,x), 
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where 8 e is the fractional energy density perturbation (we are interested in the fractional rather than 
in the absolute energy density fluctuation!), and the pressure perturbation 8p is defined as in 
In addition, there is the possibility of a non-vanishing entropy perturbation defined as in (p 



We now insert this ansatz into the basic hydrodynamical equations ( |2.2| ), linearize in the per- 
turbation variables, and combine the first order differential equations for 8 e and 8p into a single 
second order differential equation for 8p e . The result simplifies if we work in "comoving coordi- 
nates" q which are the coordinates painted onto the expanding background, i.e. x(?) = a(t)q(t) . 
After some algebra, we obtain the following equation which describes the time evolution of density 
fluctuations: 

2 

Se + 2HS e - %Vl8 e -47iG Po 8 e = -%8S, (2.10) 

where the subscript q on the V operator indicates that derivatives with respect to comoving coordi- 
nates are used. In addition, we have the equation of entropy conservation 8S = . 



Comparing with the equations (2.4) obtained in the absence of an expanding background, 
we see that the only difference is the presence of a Hubble damping term in the equation for 8 e . 
This term will moderate the exponential instability of the background to long wavelength density 
fluctuations. In addition, it will lead to a damping of the oscillating solutions on short wavelengths. 



More specifically, for physical wavenumbers k p <C kj (where kj is again given by (|2.7D), and in 



a matter-dominated background cosmology, the general solution of (|2T0J) in the absence of any 
entropy fluctuations is given by 

8 k (t) = cit 2/3 +c 2 r\ (2.11) 

where c\ and c 2 are two constants determined by the initial conditions, and we have dropped the 
subscript e in expressions involving 8 e . There are two fundamental solutions, the first is a growing 
mode with 8k(t) ~ a(t), the second a decaying mode with 8k (t) ~ t . On short wavelengths, one 
obtains damped oscillatory motion: 

8 k (t) ~ a- l l 2 {t)exp(±ic s k J dt'a-\t')) . (2.12) 

The above analysis applies for t > t eq , when we are following the fluctuations of the dominant 
component of matter. In the radiation era, cold matter fluctuations only grow logarithmically in 
time since the dominant component of matter (the relativistic radiation) does not cluster and hence 
the term in the equation of motion for fluctuations which represents the gravitational attraction 
force is suppressed. In this sense, we can say that inhomogeneities begin to cluster at the time t eq . 

As a simple application of the Newtonian equations for cosmological perturbations derived 
above, let us compare the predicted cosmic microwave background (CMB) anisotropics in a spa- 
tially flat universe with only baryonic matter - Model A - to the corresponding anisotropies in a 
flat Universe with mostly cold dark matter (pressure-less non-baryonic dark matter) - Model B. 
We start with the observationally known amplitude of the relative density fluctuations today (time 
?o), and we use the fact that the amplitude of the CMB anisotropies on the angular scale 6(k) 
corresponding to the comoving wavenumber k is set by the primordial value of the gravitational 
potential <p - introduced in the following section - which in turn is related to the primordial value 
of the density fluctuations at Hubble radius crossing (and not to its value at the time t rec - see e.g. 
Chapter 17 of [0]). 
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In Model A, the dominant component of the pressure-less matter is coupled to radiation be- 
tween t eq and t rec , the time of last scattering. Thus, the Jeans length is comparable to the Hubble 
radius. Therefore, for comoving galactic scales, k 3> kj in this time interval, and thus the frac- 
tional density contrast decreases as a(t)~ l l 2 . In contrast, in Model B, the dominant component 
of pressure-less matter couples only weakly to radiation, and hence the Jeans length is negligibly 
small. Thus, in Model B, the relative density contrast grows as a(t) between t eq and t rec . In the 
time interval t rec < t < to, the fluctuations scale identically in Models A and B. Summarizing, we 
conclude, working backwards in time from a fixed amplitude of 8k today, that the amplitudes of 
Sk{t e q) in Models A and B (and thus their primordial values) are related by 

5 k (t eq )\ A ^(^-f 2 8 k (t eq )\ B . (2.13) 
a(t eq ) 

Hence, in Model A (without non-baryonic dark matter) the CMB anisotropics are predicted to be 



a factor of about 30 larger [35 ] than in Model B, way in excess of the recent observational results. 
This is one of the strongest arguments for the existence of non-baryonic dark matter 3 . Note that 
the precise value of the enhancement factor depends on the value of the cosmological constant A - 
the above result holds for A = 0. 



2.2 Observables 

Let us consider perturbations on a fixed comoving length scale given by a comoving wavenum- 
ber k. The corresponding physical length increases as a{t). This is to be compared with the 
Hubble radius H~ l (t) which scales as t provided a(t) grows as a power of t. In the late time 
Universe, a(t) ~ f 1 / 2 in the radiation-dominated phase (i.e. for t < t eq ), and a(t) ~ ? 2 / 3 in the 
matter-dominated period {t eq < t < to). Thus, we see that at sufficiently early times, all comoving 
scales had a physical length larger than the Hubble radius. If we consider large cosmological scales 
(e.g. those corresponding to the observed CMB anisotropics or to galaxy clusters), the time tu (k) 
of "Hubble radius crossing" (when the physical length was equal to the Hubble radius) was in fact 
later than t eq . 

Cosmological fluctuations can be described either in position space or in momentum space. 
In position space, we compute the root mean square mass fluctuation 8M/M(k,i) in a sphere of 
radius / = 2%jk at time t. A scale-invariant spectrum of fluctuations is defined by the relation 

^L(k,t H (k)) = const.. (2.14) 

Such a spectrum was first suggested by Harrison [[36|] and Zeldovich [[37|] as a reasonable choice for 
the spectrum of cosmological fluctuations. We can introduce the "spectral index" n of cosmological 
fluctuations by the relation 

(^) 2 (k,t H (k))~k n -\ (2-15) 
and thus a scale-invariant spectrum corresponds to n = 1. 

3 In my opinion no proposed alternative to particle dark matter should be taken serious unless their proponents can 
demonstrate that their model can reproduce the agreement between the amplitudes of the CMB anisotropies on one hand 
and of the large-scale matter power spectrum on the other hand. 
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To make the transition to the (more frequently used) momentum space representation, we 
Fourier decompose the fractional spatial density contrast 



8 £ ( x ,t) =V l/2 J d 3 k8 e (k,t)e ikx , 



where V is a cutoff volume. The power spectrum P§ of density fluctuations is defined by 

P s (k) =k 3 \5 e (k)\ 2 , 



(2.16) 



(2.17) 



where k is the magnitude of k, and we have assumed for simplicity a Gaussian distribution of 
fluctuations in which the amplitude of the fluctuations only depends on k. 
We can also define the power spectrum of the gravitational potential (p: 



P 9 (k) = k 3 \8(p(k)\ 2 . 



These two power spectra are related by the Poisson equation ( |2.2| ) 



P v (k) ~ k- 4 P s (k). 



(2.18) 



(2.19) 



In general, the condition of scale-invariance is expressed in momentum space in terms of the 
power spectrum evaluated at a fixed time. To obtain this condition, we first use the time dependence 



of the fractional density fluctuation from ( [2.1 1[ ) to determine the mass fluctuations at a fixed time 
t > tji(k) > t eq (the last inequality is a condition on the scales considered) 



< *-) (M) = ( ^) 



.4/3 



,8M, 



~(k,t H (k)). 



The time of Hubble radius crossing is given by 

a{t H {k))k- 1 = 2t H {k), 

and thus 



t H (k) 



1/3 



Inserting this result into ( |2.20| ) and making use of ( [2.1 5\ > we find 



(^) 2 (M)~*"«. 



(2.20) 



(2.21) 



(2.22) 



(2.23) 



Since, for reasonable values of the index of the power spectrum, 8M/M(k,t) is dominated by the 

(2.24) 



Fourier modes with wavenumber k, we find that (2.23) implies 



|4| 2 ~*' ! , 



or, equivalently, 



P 9 (k) ~ k 



n-\ 



(2.25) 
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2.3 Classical Relativistic Theory 
2.3.1 Introduction 

The Newtonian theory of cosmological fluctuations discussed in the previous section breaks 
down on scales larger than the Hubble radius because it neglects perturbations of the metric, and 
because on large scales the metric fluctuations dominate the dynamics. 

Let us begin with a heuristic argument to show why metric fluctuations are important on scales 
larger than the Hubble radius. For such inhomogeneities, one should be able to approximately 



describe the evolution of the space-time by applying the first FRW equation (3.7) of homogeneous 
and isotropic cosmology to the local Universe (this approximation is made more rigorous in [38]). 
Based on this equation, a large-scale fluctuation of the energy density will lead to a fluctuation 
("8a") of the scale factor a which grows in time. This is due to the fact that self gravity amplifies 
fluctuations even on length scales A greater than the Hubble radius. 

This argument is made rigorous in the following analysis of cosmological fluctuations in the 
context of general relativity, where both metric and matter inhomogeneities are taken into account. 
We will consider fluctuations about a homogeneous and isotropic background cosmology, given by 



the metric (IT), which can be written in conformal time r\ (defined by dt = a(t)dr\) as 

ds 2 = a(r)) 2 {dr} 2 -dx 2 ). (2.26) 

The theory of cosmological perturbations is based on expanding the Einstein equations to lin- 
ear order about the background metric. The theory was initially developed in pioneering works 



by Lifshitz [39]. Significant progress in the understanding of the physics of cosmological fluctu- 



ations was achieved by Bardeen [ 4-0 ] who realized the importance of subtracting gauge artifacts 



(see below) from the analysis (see also [41, |2|]). The following discussion is based on Part I of the 
comprehensive review article [|4j]. Other reviews - in some cases emphasizing different approaches 
-are @,00|§]. 

2.3.2 Classifying Fluctuations 

The first step in the analysis of metric fluctuations is to classify them according to their trans- 
formation properties under spatial rotations. There are sc alar, vector and second rank tensor fluc- 
tuations. In linear theory, there is no coupling between the different fluctuation modes, and hence 
they evolve independently (for some subtleties in this classification, see [|47|]). 



We begin by expanding the metric about the FRW background metric g^l given by ( |2.26| ): 



guv = sfl + 8g^ v ■ (2.27) 

The background metric depends only on time, whereas the metric fluctuations 8g^ v depend on 
both space and time. Since the metric is a symmetric tensor, there are at first sight 10 fluctuating 
degrees of freedom in 8g^ v . 

There are four degrees of freedom which correspond to scalar metric fluctuations (the only 
four ways of constructing a metric from scalar functions): 



18 



Cosmology 



Robert H. Brandenberger 



where the four fluctuating degrees of freedom are denoted (following the notation of [Q]) <p,B,E, 
and a comma denotes the ordinary partial derivative (if we had included spatial curvature of the 
background metric, it would have been the covariant derivative with respect to the spatial metric), 
and Sij is the Kronecker symbol. 

There are also four vector degrees of freedom of metric fluctuations, consisting of the four 
ways of constructing metric fluctuations from three vectors: 




5 *f v =a Z \ c r ' , (2-29) 



where S; and Fj are two divergence-less vectors (for a vector with non- vanishing divergence, the 
divergence contributes to the scalar gravitational fluctuation modes). 

Finally, there are two tensor modes which correspond to the two polarization states of gravita- 
tional waves: 

S 8liV = -a 2 I ° ,° J , (2.30) 



hi 



where hn is trace-free and divergence-less 



h\ = hjj = 0. (2.31) 

Gravitational waves do not couple at linear order to the matter fluctuations. Vector fluctu- 
ations decay in an expanding background cosmology and hence are not usually cosmologically 
important. The most important fluctuations, at least in inflationary cosmology, are the scalar metric 
fluctuations, the fluctuations which couple to matter inhomogeneities and which are the relativistic 
generalization of the Newtonian perturbations considered in the previous section. Note that vector 



and tensor perturbations are important (see e.g. [|48|]) in topological defects models of structure 



formation and in bouncing cosmologies [19] 



2.3.3 Gauge Transformation 

The theory of cosmological perturbations is at first sight complicated by the issue of gauge 
invariance (at the final stage, however, we will see that we can make use of the gauge freedom to 
substantially simplify the theory). The coordinates t,x of space-time carry no independent physical 
meaning. They are just labels to designate points in the space-time manifold. By performing a 
small-amplitude transformation of the space-time coordinates (called "gauge transformation" in 
the following), we can easily introduce "fictitious" fluctuations in a homogeneous and isotropic 
Universe. These modes are "gauge artifacts". 

We will in the following take an "active" view of gauge transformation. Let us consider two 
space-time manifolds, one of them a homogeneous and isotropic Universe ^o, the other a physical 
Universe with inhomogeneities. A choice of coordinates can be considered to be a mapping 9l 
between the manifolds and ^ . Let us consider a second mapping Qi which will take the same 
point (e.g. the origin of a fixed coordinate system) in to a different point in Using the 
inverse of these maps and we can assign two different sets of coordinates to points in jjt, '. 

Consider now a physical quantity Q (e.g. the Ricci scalar) on and the corresponding 
physical quantity Q^> on Then, in the first coordinate system given by the mapping the 
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perturbation 8Q of Q at the point p £ ^ is defined by 

52b) = G(P)-G (0) (^ _1 (P))- ( 2 -32) 
Analogously, in the second coordinate system given by Q>, the perturbation is defined by 

5Q(p) = Q(p)-Q {0) (@- 1 (p)). (2.33) 

The difference 

Afi(p) = S~<2(/>) - 5fi(p) (2.34) 

is obviously a gauge artifact and carries no physical significance. 

Some of the degrees of freedom of metric perturbation introduced in the first subsection are 
gauge artifacts. To isolate these, we must study how coordinate transformations act on the metric. 
There are four independent gauge degrees of freedom corresponding to the coordinate transforma- 
tion 

+ (2.35) 

The zero (time) component £° of ^ leads to a scalar metric fluctuation. The spatial three vector 
E, 1 can be decomposed as 

V & ■ Y'li (2-36) 

(where is the spatial background metric) into a transverse piece Q r which has two degrees of 
freedom which yield vector perturbations, and the second term (given by the gradient of a scalar 
which gives a scalar fluctuation. To summarize this paragraph, there are two scalar gauge modes 
given by £° and and two vector modes given by the transverse three vector Q r . Thus, there 
remain two physical scalar and two vector fluctuation modes. The gravitational waves are gauge- 
invariant. 

Let us now focus on how the scalar gauge transformations (i.e. the transformations given by 
£° and act on the scalar metric fluctuation variables <j),B,E, and \ff. An immediate calculation 
yields: 

B = B + ^-t,' (2.37) 

E = E-% 

a' n 
a 

where a prime indicates the derivative with respect to conformal time X] . 

There are two approaches to deal with the gauge ambiguities. The first is to fix a gauge, 
i.e. to pick conditions on the coordinates which completely eliminate the gauge freedom, the 
second is to work with a basis of gauge-invariant variables. If one wants to adopt the gauge-fixed 
approach, there are many different gauge choices. Note that the often used synchronous gauge 
determined by <5g° M = does not totally fix the gauge. A convenient system which completely 
fixes the coordinates is the so-called longitudinal or conformal Newtonian gauge defined by 
B = E = 0. This is the gauge we will use in the following. 
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Note that the above analysis was strictly at linear order in perturbation theory. Beyond linear 
order, the structure of perturbation theory becomes much more involved. In fact, one can show 
[ |50| ] that the only fluctuation variables which are invariant under all coordinate transformations are 
perturbations of variables which are constant in the background space-time. 

2.3.4 Equations of Motion 

We begin with the Einstein equations 

G MV = 87tGT^ v , (2.38) 

where G^ v is the Einstein tensor associated with the space-time metric g^ v , and 7^ v is the energy- 
momentum tensor of matter, insert the ansatz for metric and matter perturbed about a FRW back- 
ground (gjS(Tj), <P (0) C7)): 

S/xvM) = gflW + Sg^r,) (2.39) 
<p(x,T}) = <po(T?) + <5<p(x,T7), (2.40) 

(where we have for simplicity replaced general matter by a scalar matter field (p) and expand to 
linear order in the fluctuating fields, obtaining the following equations: 

<5G MV = SxGdT^ . (2.41) 

In the above, dg^ v is the perturbation in the metric and S(p is the fluctuation of the matter field (p. 
Inserting the longitudinal gauge metric 

ds 2 = a 2 [(1 + 2§)dr} 2 - (1 - ly^ijdjdxi] (2.42) 

into the general perturbation equations ( [2.4 1[ ) yields the following set of equations of motion: 

-3jr(j^0 + y/) +VV = 4xGa 2 ST$ 
(jSTQ + yr) i = 4xGa 2 8T? 
[[2je' + J^ 2 )(j) + Jf<p' + y" + 2^1//] 8) (2.43) 

+ l -^DS)- l -Y k D, kj = —AnGa 2 STj , 

where D = (j) — y/ and M J = a' /a. 

The first conclusion we can draw is that if no anisotropic stress is present in the matter at linear 
order in fluctuating fields, i.e. 8T] = for i ^ j, then the two metric fluctuation variables coincide: 

= i//\ (2.44) 

This will be the case in most simple cosmological models, e.g. in theories with matter described 
by a set of scalar fields with canonical form of the action, and in the case of a perfect fluid with no 
anisotropic stress. 

Let us now restrict our attention to the case of matter described in terms of a single scalar 
field cp (in which case the perturbations on large scales are automatically adiabatic) which can be 
expanded as 

9(x,T]) = <p (T]) + 5<p(x,T]) (2.45) 
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in terms of background matter <po and matter fluctuation 8(p(x,rj). Then, in longitudinal gauge, 



( j2-43[ ) reduce to the following set of equations of motion 



V 2 </>-3jr</>'-(jr'+2^ 2 )0 = 4nG(%8(p' + V'a 2 8(p) 
</>' + = 47iG(p' 8(p 
0" + 3^0' + (jf'+2^ 2 )0 = 4nG(%Sq>' -v'a 2 Sq>) 



(2.46) 



where V denotes the derivative of V with respect to (p. These equations can be combined to give 
the following second order differential equation for the relativistic potential <p : 



(2.47) 



This is the final result for the classical evolution of cosmological fluctuations. First of all, we 



note the similarities with the equation (2.10) obtained in the Newtonian theory. The final term in 



( |2.47[ ) is the force due to gravity leading to the instability, the second to last term is the pressure 
force leading to oscillations (relativistic since we are considering matter to be a relativistic field), 
and the second term represents the Hubble friction. For each wavenumber there are two funda- 
mental solutions. On small scales (k > H), the solutions correspond to damped oscillations, on 
large scales (k < H) the oscillations freeze out and the dynamics is governed by the gravitational 
force competing with the Hubble friction term. Note, in particular, how the Hubble radius naturally 
emerges as the scale where the nature of the fluctuating modes changes from oscillatory to frozen. 

Considering the equation in a bit more detail, we observe that if the equation of state of the 
background is independent of time (which will be the case if =%= 0), then in an expanding 



background, the dominant mode of (2.47) is constant, and the sub-dominant mode decays. If the 
equation of state is not constant, then the dominant mode is not constant in time. Specifically, at the 
end of inflation Jf? < 0, and this leads to a growth of (see the following subsection). In contrast, 
in a contracting phase the dominant mode of (p is growing. 



To study the quantitative implications of the equation of motion ( |2.47| ), it is convenient to 
introduce [51 , |2|| the variable £ (which, up to a correction term of the order V 2 </> which is unim- 
portant for large-scale fluctuations, is equal to the curvature perturbation in comoving gauge 



[53]) by 



2 (g- 1 J + ») 
V 3 l+w 



where 



w = — 



(2.48) 



(2.49) 



characterizes the equation of state of matter. In terms of £, the equation of motion ( |2.47| ) takes on 
the form 



<ff(V 2 <$>). 



(2.50) 



On large scales, the right hand side of the equation is negligible, which leads to the conclusion that 
large-scale cosmological fluctuations satisfy 



f(l + w) = 0. 



(2.51) 
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This implies that £ is constant except possibly if 1 + w = at some point in time during the 
cosmological evolution (which occurs during reheating in inflationary cosmology if the inflaton 



field undergoes oscillations [|54|]). In single matter field models it is indeed possible to show that 
£ = on super-Hubble scales independent of assumptions on the equation of state [ j55j |56| ]. This 
"conservation law" makes it easy to relate initial fluctuations to final fluctuations in inflationary 
cosmology, as will be illustrated in the following subsection. 

In the presence of non-adiabatic fluctuations (entropy fluctuations) there is a source term on the 



righ-hand side of (2.51) which is proportional to the amplitude of the entropy fluctuations. Thus, 
as already seen in the case of the Newtonian theory of cosmological fluctuations, non-adiabatic 
fluctuations induce a growing adiabatic mode (see [ 57[ 58] for discussions of the consequences in 
double field inflationary models). 

2.4 Quantum Theory 
2.4.1 Overview 

As we have seen in the first section, in many models of the very early Universe, in particular 
in inflationary cosmology, in string gas cosmology and in the matter bounce scenario, primordial 
inhomogeneities are generated in an initial phase on sub-Hubble scales. The wavelength is then 
stretched relative to the Hubble radius, becomes larger than the Hubble radius at some time and 
then propagates on super-Hubble scales until re-entering at late cosmological times. In a majority 
of the current structure formation scenarios (string gas cosmology is an exception in this respect), 
fluctuations are assumed to emerge as quantum vacuum perturbations. Hence, to describe the 
generation and evolution of the inhomogeneities a quantum treatment is required. 

In the context of a Universe with a de Sitter phase, the quantum origin of cosmological fluc- 
tuations was first discussed in ^ and also [J7| for earlier ideas. In particular, Mukhanov [||] 
and Press [|7|] realized that in an exponentially expanding background, the curvature fluctuations 
would be scale-invariant, and Mukhanov provided a quantitative calculation which also yielded the 
logarithmic deviation from exact scale-invariance. 

The role of the Hubble radius has already been mentioned repeatedly in these lectures. In 
particular, in the previous subsection we saw that the Hubble radius separates scales on which fluc- 
tuations oscillate (sub-Hubble scales) from those where they are frozen in (super-Hubble scales). 
Another way to see the role of the Hubble radius is to consider the equation of a free scalar matter 
field q> on an unperturbed expanding background: 

V 2 

(p + 3Hw-^-(p = 0. (2.52) 

or 

The second term on the left hand side of this equation leads to damping of cp with a characteristic 
decay rate given by H. As a consequence, in the absence of the spatial gradient term, (p would be 
of the order of magnitude Hq>. Thus, comparing the second and the third terms on the left hand 
side, we immediately see that the microscopic (spatial gradient) term dominates on length scales 
smaller than the Hubble radius, leading to oscillatory motion, whereas this term is negligible on 
scales larger than the Hubble radius, and the evolution of cp is determined primarily by gravity. 
Note that in general cosmological models the Hubble radius is much smaller than the horizon (the 
forward light cone calculated from the initial time). In an inflationary universe, the horizon is 
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larger by a factor of at least exp(Af), where N is the number of e-foldings of inflation, and the 
lower bound is taken on if the Hubble radius and horizon coincide until inflation begins. It is very 
important to realize this difference, a difference which is obscured in most articles on cosmology 
in which the term "horizon" is used when "Hubble radius" is meant. Note, in particular, that the 
homogeneous inflaton field contains causal information on super-Hubble but sub-horizon scales. 
Hence, it is completely consistent with causality [|54j] to have a microphysical process related to the 
background scalar matter field lead to exponential amplification of the amplitude of fluctuations 
during reheating on such scales, as it does in models in which entropy perturbations are present 



and not suppressed during inflation [ ]57| , p8| ]. Note that also in string gas cosmology and in the 
matter bounce scenario the Hubble radius and horizon are completely different. 

There are general relativistic conservation laws [[>9|] which imply that adiabatic fluctuations 
produced locally must be Poisson-statistic suppressed on scales larger than the Hubble radius. For 
example, fluctuations produced by the formation of topological defects at a phase transition in the 



early universe are initially isocurvature (entropy) in nature (see e.g. [50] for a discussion). Via the 



source term in the equation of motion (2.4), a growing adiabatic mode is induced, but at any fixed 



time the spectrum of the curvature fluctuation on scales larger than the Hubble radius has index 



n = 4 (Poisson). A similar conclusion applies to models [ pl| , p2| | of modulated reheating (see [ p3[ ] 
for a nice discussion), and to models in which moduli fields obtain masses after some symmetry 
breaking, their quantum fluctuations then inducing cosmological perturbations. A prototypical 
example is given by axion fluctuations in an inflationary universe (see e.g. [ ]64| ] and references 
therein). 

To understand the generation and evolution of fluctuations in current models of the very early 
Universe, we need both Quantum Mechanics and General Relativity, i.e. quantum gravity. At first 
sight, we are thus faced with an intractable problem, since the theory of quantum gravity is not 
yet established. We are saved by the fact that today on large cosmological scales the fractional 
amplitude of the fluctuations is smaller than 1. Since gravity is a purely attractive force, the fluctu- 
ations had to have been - at least in the context of an eternally expanding background cosmology 
- very small in the early Universe. Thus, a linearized analysis of the fluctuations (about a classical 
cosmological background) is self-consistent. 

From the classical theory of cosmological perturbations discussed in the previous section, we 
know that the analysis of scalar metric inhomogeneities can be reduced - after extracting gauge 
artifacts - to the study of the evolution of a single fluctuating variable. Thus, we conclude that the 
quantum theory of cosmological perturbations must be reducible to the quantum theory of a single 
free scalar field which we will denote by v. Since the background in which this scalar field evolves 
is time-dependent, the mass of v will be time-dependent. The time-dependence of the mass will 
lead to quantum particle production over time if we start the evolution in the vacuum state for v. 
As we will see, this quantum particle production corresponds to the development and growth of 
the cosmological fluctuations. Thus, the quantum theory of cosmological fluctuations provides a 
consistent framework to study both the generation and the evolution of metric perturbations. The 



following analysis is based on Part II of [ 34 ] . 

2.4.2 Outline of the Analysis 

In order to obtain the action for linearized cosmological perturbations, we expand the action 
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to quadratic order in the fluctuating degrees of freedom. The linear terms cancel because the 
background is taken to satisfy the background equations of motion. 

We begin with the Einstein-Hilbert action for gravity and the action of a scalar matter field (for 



the more complicated case of general hydrodynamical fluctuations the reader is referred to [p4||) 



S = Jd*xV=g[-j^R + ^d lt <pd' l <p-V(<p)] 



(2.53) 



l6nG 2 

where R is the Ricci curvature scalar. 

The simplest way to proceed is to work in longitudinal gauge, in which the metric and matter 
take the form (assuming no anisotropic stress) 



ds 2 = a 2 (T7)[(l+20(T7,x))dT7 2 
9(77, x) = (po(T]) + 5<p(77,x). 



(1 -2<l>(t,x))dx 2 



(2.54) 



The two fluctuation variables and (p must be linked by the Einstein constraint equations since 
there cannot be matter fluctuations without induced metric fluctuations. 



The two nontrivial tasks of the lengthy Q34[ ] computation of the quadratic piece of the action 
is to find out what combination of (p and gives the variable v in terms of which the action 
has canonical kinetic term, and what the form of the time-dependent mass is. This calculation 



involves inserting the ansatz ( |2.54[ ) into the action ( |2.53[ ), expanding the result to second order in 
the fluctuating fields, making use of the background and of the constraint equations, and dropping 
total derivative terms from the action. In the context of scalar field matter, the quantum theory of 
cosmological fluctuations was developed by Mukhanov [^, |66|] and Sasaki [67]. The result is the 
following contribution to the action quadratic in the perturbations: 



5(2) 



d 4 x[v 



J2 



Z" 21 

•v/v.H v 



(2.55) 

where the canonical variable v (the "Sasaki-Mukhanov variable" introduced in [Q - see also [|68|]) 
is given by 



a[Sq> + ^], 



with ,%f = d ja, and where 



a% 



(2.56) 



(2.57) 



As long as the equation of state does not change over time) 

Note that the variable v is related to the curvature perturbation «■ 
duced in [ 53[ | and closely related to the variable £ used in & 

v = 



(2.58) 

in comoving coordinates intro- 



(2.59) 



The equation of motion which follows from the action ( |2.55| ) is (in momentum space) 



"2 ^ 
v k +k v k v k = 0, 



(2.60) 
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where v& is the k'th Fourier mode of v. The mass term in the above equation is in general given 
by the Hubble scale. Thus, it immediately follows that on small length scales, i.e. for k > Ich, the 
solutions for Vk are constant amplitude oscillations . These oscillations freeze out at Hubble radius 
crossing, i.e. when k = ku- On longer scales (k <C there is a mode of which scales as z. This 
mode is the dominant one in an expanding universe, but not in a contracting one. 

Given the action ( 2.55 ), the quantization of the cosmological perturbations can be performed 
by canonical quantization (in the same way that a scalar matter field on a fixed cosmological back- 
ground is quantized [|6^]). 

The final step in the quantum theory of cosmological perturbations is to specify an initial 
state. Since in inflationary cosmology all pre-existing classical fluctuations are red-shifted by the 
accelerated expansion of space, one usually assumes (we will return to a criticism of this point 
when discussing the trans-Planckian problem of inflationary cosmology) that the field v starts out 
at the initial time t{ mode by mode in its vacuum state. Two questions immediately emerge: what 
is the initial time t{, and which of the many possible vacuum states should be chosen. It is usually 
assumed that since the fluctuations only oscillate on sub-Hubble scales, the choice of the initial 
time is not important, as long as it is earlier than the time when scales of cosmological interest 
today cross the Hubble radius during the inflationary phase. The state is usually taken to be the 
Bunch-Davies vacuum (see e.g. [^D, since this state is empty of particles at t\ in the coordinate 
frame determined by the FRW coordinates Thus, we choose the initial conditions 

v k {m) = -L (2.61) 
\/2k 



, \[k 

where tj, is the conformal time corresponding to the physical time 
Returning to the case of an expanding universe, the scaling 

vt ~ z ~ a (2.62) 

implies that the wave function of the quantum variable Vk which performs quantum vacuum fluc- 
tuations on sub-Hubble scales, stops oscillating on super-Hubble scales and instead is squeezed 
(the amplitude increases in configuration space but decreases in momentum space). This squeezing 
corresponds to quantum particle production. It is also one of the two conditions which are required 
for the classicalization of the fluctuations. The second condition is decoherence which is induced 
by the non-linearities in the dynamical system which are inevitable since the Einstein action leads 
to highly nonlinear equatiions (see [f7(J for a recent discussion of this point, and [ f7l| ] for related 
work). 

Note that the squeezing of cosmological fluctuations on super-Hubble scales occurs in all 
models, in particular in string gas cosmology and in the bouncing universe scenario since also in 
these scenarios perturbations propagate on super-Hubble scales for a long period of time. In a 
contracting phase, the dominant mode of on super-Hubble scales is not the one given in ( 2.62 ) 
(which in this case is a decaying mode), but the second mode which scales as z~ p with an exponent 
p which is positive and whose exact value depends on the background equation of state. 

Applications of this theory in inflationary cosmology, in the matter bounce scenario and in 
string gas cosmology will be considered in the respective sections of these lecture notes. 
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2.4.3 Quantum Theory of Gravitational Waves 

The quantization of gravitational waves parallels the quantization of scalar metric fluctuations, 
but is more simple because there are no gauge ambiguities. Note that at the level of linear fluctua- 
tions, scalar metric fluctuations and gravitational waves are independent. Both can be quantized on 
the same cosmological background determined by the background scale factor and the background 
matter. However, in contrast to the case of scalar metric fluctuations, the tensor modes are also 
present in pure gravity (i.e. in the absence of matter). 

Starting point is the action ( [2.53 ). Into this action we insert the metric which corresponds to a 
classical cosmological background plus tensor metric fluctuations: 

ds 2 = a 2 (ri)[dri 2 - (<% + h^dx'd^] , (2.63) 

where the second rank tensor fcy(7],x) represents the gravitational waves, and in turn can be de- 
composed as 

hij(T],x) = h + (Ti,x)efj+h x (Ti,x)eXj (2.64) 

into the two polarization states. Here, e^- and efj are two fixed polarization tensors, and h + and h x 
are the two coefficient functions. 

To quadratic order in the fluctuating fields, the action consists of separate terms involving h + 
and h x . Each term is of the form 

2 

S W = J d 4 x^-[ti 2 - (Vh) 2 ] , (2.65) 

leading to the equation of motion 

h" k + + k 2 h k = . (2.66) 



The variable in terms of which the action (2.65) has canonical kinetic term is 

H = ah k , (2.67) 

and its equation of motion is 

A4' + (* 2 -^)At* = 0. (2.68) 



This equation is very similar to the corresponding equation (2.60) for scalar gravitational inhomo- 
geneities, except that in the mass term the scale factor a(r\) replaces z(tj), which leads to a very 
different evolution of scalar and tensor modes during the reheating phase in inflationary cosmology 
during which the equation of state of the background matter changes dramatically. 

Based on the above discussion we have the following theory for the generation and evolution 
of gravitational waves in an accelerating Universe (first developed by Grishchuk [|72|]): waves exist 
as quantum vacuum fluctuations at the initial time on all scales. They oscillate until the length 
scale crosses the Hubble radius. At that point, the oscillations freeze out and the quantum state of 
gravitational waves begins to be squeezed in the sense that 

ttfc(T?)~fl(n), (2-69) 



which, from (2.67) corresponds to constant amplitude of h k . The squeezing of the vacuum state 
leads to the emergence of classical properties of this state, as in the case of scalar metric fluctua- 
tions. 
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Figure 9: Space-time diagram (sketch) showing the evolution of scales in inflationary cosmology. The 
vertical axis is time, and the period of inflation lasts between f; and tg, and is followed by the radiation- 
dominated phase of standard big bang cosmology. During exponential inflation, the Hubble radius H^ 1 is 
constant in physical spatial coordinates (the horizontal axis), whereas it increases linearly in time after 
The physical length corresponding to a fixed comoving length scale labelled by its wavenumber k increases 
exponentially during inflation but increases less fast than the Hubble radius (namely as f 1 ' 2 ), after inflation. 

3. Inflationary Cosmology 

3.1 Mechanism of Inflation 

The idea of inflationary cosmology is to assume that there was a period in the very early 
Universe during which the scale factor was accelerating, i.e. a > 0. This implies that the Hubble 
radius was shrinking in comoving coordinates, or, equivalently, that fixed comoving scales were 
"exiting" the Hubble radius. In the simplest models of inflation, the scale factor increases nearly 
exponentially. As illustrated in Figure ^j, the basic geometry of inflationary cosmology provides 
a solution of the fluctuation problem. As long as the phase of inflation is sufficiently long, all 
length scales within our present Hubble radius today originate at the beginning of inflation with a 
wavelength smaller than the Hubble radius at that time. Thus, it is possible to create perturbations 
locally using physics obeying the laws of special relativity (in particular causality). As will be 
discussed later, it is quantum vacuum fluctuations of matter fields and their associated curvature 
perturbations which are responsible for the structure we observe today. 

Postulating a phase of inflation in the very early universe solves the horizon problem of the 
SBB, namely it explains why the causal horizon at the time t rec when photons last scatter is larger 
than the radius of the past light cone at t rec , the part of the last scattering surface which is visible 
today in CMB experiments. Inflation also explains the near flatness of the universe: in a decelerat- 
ing universe spatial flatness is an unstable fixed point of the dynamics, whereas in an accelerating 
universe it becomes an attractor. Another important aspect of the inflationary solution of the flat- 
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ness problem is that inflation exponentially increases the volume of space. Without inflation, it is 
not possible that a Planck scale universe at the Planck time evolves into a sufficiently large universe 
today. 

3.2 Fluctuations in Inflationary Cosmology 

We will now use the quantum theory of cosmological perturbations developed in the previous 
section to calculate the spectrum of curvature fluctuations in inflationary cosmology. The starting 
point are quantum vacuum initial conditions for the canonical fluctuation variable v k : 

1 



2k 



(3.1) 



for all k for which the wavelength is smaller than the Hubble radius at the initial time t{. 

The amplitude remains unchanged until the modes exit the Hubble radius at the respective 
times tn (k) given by 



a-\t H {k))k = H. 

We need to compute the power spectrum £Pgg(k) of the curvature fluctuation 



(3.2) 
defined in 



( |2.59| ) at some late time t when the modes are super-Hubble. We first relate the power spectrum via 



the growth rate ( 2.62 ) of v on super-Hubble scales to the power spectrum at the time tn{k) and then 
use the constancy of the amplitude of v on sub-Hubble scales to relate it to the initial conditions 
©. Thus 

k^ 2 k (t) = k\- 2 {t)\v k {t)\ 2 



f (k,t) 



(3.3) 



k\-\t){ 



a(t) 



k"z 



a{t H {k)) 
2 {t H (k))\v k (t H {k))\ 2 



v k (t H (k))\ : 



k\^) 2 a- 2 (t H (k))\v k ( ti )\ 2 . 



where in the final step we have used (2.58) and the constancy of the amplitude of v on sub-Hubble 



scales. Making use of the condition (3.2) for Hubble radius crossing, and of the initial conditions 



(3.1), we immediately see that 



' a{t) ^k- 2 k-'H 2 



z(t) 



(3.4) 



and that thus a scale invariant power spectrum with amplitude proportional to H 2 results, in agree- 
ment with what was argued on heuristic grounds in the overview of inflation in the the first section. 
To obtain the precise amplitude, we need to make use of the relation between z and a. We obtain 

H 4 



(3.5) 



which for any given value of k is to be evaluated at the time tn{k) (before the end of inflation). For 
a scalar field potential (see following subsection) 

v 



V(q>) 



(3.6) 



the resulting amplitude in (3.5) is X. Thus, in order to obtain the observed value of the power 
spectrum of the order of 10~ 10 , the coupling constant A must be tuned to a very small value. 
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3.3 Models of Inflation 

Let us now consider how it is possible to obtain a phase of cosmological inflation. We will 
assume that space-time is described using the equations of General Relativity 4 . In this case, 
the dynamics of the scale factor a(t) is determined by the Friedmann-Robertson- Walker (FRW) 
equations 

(-) 2 = 87lGp (3.7) 
a 

and 

t = -*Jf {p + 3 p) (3.8) 

where for simplicity we have omitted the contributions of spatial curvature (since spatial curvature 
is diluted during inflation) and of the cosmological constant (since any small cosmological con- 
stant which might be present today has no effect in the early Universe since the associated energy 
density does not increase when going into the past). From (^8) it is clear that in order to obtain an 
accelerating universe, matter with sufficiently negative pressure 

P < ~\p (3-9) 

is required. Exponential inflation is obtained for p = —p. 

Conventional perfect fluids have positive semi-definite pressure and thus cannot yield inflation. 
Quantum field theory can come to the rescue. We know that a description of matter in terms of 
classical perfect fluids must break down at early times. An improved description of matter will be 
given in terms of quantum fields. Matter which we observe today consists of spin 1/2 and spin 1 
fields. Such fields cannot yield inflation in the context of the renormalizable quantum field theory 
framework. The existence of scalar matter fields is postulated to explain the generation of mass via 
spontaneous symmetry breaking. Scalar matter fields (denoted by (p) are special in that they allow 
- at the level of a renormalizable action - the presence of a potential energy term V((p). The energy 
density and pressure of a scalar field q> with canonically normalized action 5 

% = V^gi-d^d^- V(q>)] (3.10) 
(where Greek indices are space-time indices and g is the determinant of the metric) are given by 

p = l -{<p) 2 + l -a- 2 {V<p) 2 + V(<p) 

P = \(<P) 2 -l a ~ 2 (V<P) 2 -V(<p)- (3.11) 

Thus, it is possible to obtain an almost exponentially expanding universe provided the scalar field 
configuration 6 satisfies 

i(V p <p) 2 « V(q>), (3.12) 

i(<p) 2 « V(<p). (3.13) 

4 Note, however, that the first model of exponential expansion of space [Q] made use of a higher derivative gravita- 
tional action. 

5 See [j73[[74|] for discussions of fields with non-canonical kinetic terms. 
6 The scalar field yielding inflation is called the inflaton. 
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In the above, V„ = a~ l V is the gradient with respect to physical as opposed to comoving coordi- 
nates. Since spatial gradients redshift as the universe expands, the first condition will (for single 
scalar field models) always be satisfied if it is satisfied at the initial time 1 . It is the second condition 
which is harder to satisfy In particular, this condition is in general not preserved in time even it is 
initially satisfied. 

In his original model [@], Guth assumed that inflation was generated by a scalar field (p sitting 
in a false vacuum with positive vacuum energy. Hence, the conditions ( |3.12| ) and ( 3.13 ) are au- 



tomatically satisfied. Inflation ends when <p tunnels to its true vacuum state with (by assumption) 
vanishing potential energy. However, as Guth realized himself, this model is plagued by a "graceful 
exit" problem: the nucleation produces a bubble of true vacuum whose initial size is microscopic 



[J76p (see also [ |77| ] for a review of field theory methods which are useful in inflationary cosmol- 
ogy). After inflation, this bubble will expand at the speed of light, but it is too small to explain the 
observed size of the universe. Since space between the bubbles expands exponentially, the proba- 
bility of bubble percolation is negligible, at least in Einstein gravity. For this reason, the focus in 
inflationary model building shifted to the "slow-roll" paradigm [|7q1. 

It is sufficient to obtain a period of cosmological inflation that the slow-roll conditions for q> 
are satisfied. Recall that the equation of motion for a homogeneous scalar field in a cosmological 



space-time is (as follows from ( |3.10| )) is 

<p + 3H(p = -V'((p), (3.14) 

where a prime indicates the derivative with respect to (p. In order that the scalar field roll slowly, it 
is necessary that 

ip < 3Hq> (3.15) 



such that the first term in the scalar field equation of motion ( 3.14 ) is negligible. In this case, the 



condition (3.13) becomes 



V' 

(— f < 4StcG (3.16) 



and ( 3.15 ) becomes 

V" 

— « IAtzG. (3.17) 

There are many models of scalar field-driven slow-roll inflation. Many of them can be divided 
into three groups: small-field inflation, large-field inflation and hybrid inflation. Small-field infla- 
tionary models are based on ideas from spontaneous symmetry breaking in particle physics. We 
take the scalar field to have a potential of the form 

V(«p) = l -X{<p 2 -o 2 )\ (3.18) 

where a can be interpreted as a symmetry breaking scale, and A is a dimensionless coupling con- 
stant. The hope of initial small-field models ("new inflation" [f78|]) was that the scalar field would 
begin rolling close to its symmetric point <p = 0, where thermal equilibrium initial conditions would 



7 In fact, careful studies [^] show that since the gradients decrease even in a non-inflationary backgrounds, they can 
become subdominant even if they are not initially subdominant. 
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localize it in the early universe. At sufficiently high temperatures, <p = is a stable ground state of 
the one-loop finite temperature effective potential Vt(<P) (see e.g. [ |77| ] for a review). Once the tem- 
perature drops to a value smaller than the critical temperature T c , (p = turns into an unstable local 
maximum of VV(<p), and cp is free to roll towards a ground state of the zero temperature potential 
( [3.18 ). The direction of the initial rolling is triggered by quantum fluctuations. The reader can eas- 
ily check that for the potential ( [3.18 ) the slow-roll conditions cannot be satisfied if a <C m p i, where 
nipt is the Planck mass which is related to G. If the potential is modified to a Coleman- Weinberg 
@ form 

V(<p) = ^>M_l ]+ _l Av 4 ( 3. 19 ) 

(where v denotes the value of the minimum of the potential) then the slow-roll conditions can be 
satisfied. However, this corresponds to a severe fine-tuning of the shape of the potential. A further 
problem for most small-field models of inflation (see e.g. ffifo for a review) is that in order to 
end up close to the slow-roll trajectory, the initial field velocity must be constrained to be very 
small. This initial condition problem of small-field models of inflation effects a number of recently 
proposed brane inflation scenarios, see e.g. [[H]] for a discussion. 

There is another reason for abandoning small-field inflation models: in order to obtain a suffi- 
ciently small amplitude of density fluctuations, the interaction coefficients of (p must be very small. 
This makes it inconsistent to assume that <p started out in thermal equilibrium [p^]. In the absence 
of thermal equilibrium, the phase space of initial conditions is much larger for large values of (p. 

This brings us to the discussion of large-field inflation models, initially proposed in [83] under 
the name "chaotic inflation". The simplest example is provided by a massive scalar field with 
potential 

V(q>) = \m 2 <p 2 , (3.20) 

where m is the mass. It is assumed that the scalar field rolls towards the origin from large values 
of |<p|. It is a simple exercise for the reader to verify that the slow-roll conditions (3.16) and (3.17) 
are satisfied provided 

\9\>^m pl . (3.21) 

Values of | <p | comparable or greater than m p \ are also required in other realizations of large-field 
inflation. Hence, one may worry whether such a toy model can consistently be embedded in a re- 
alistic particle physics model, e.g. supergravity. In many such models V((p) receives supergravity- 
induced correction terms which destroy the flatness of the potential for \(p\ > m p i. As can be seen 
by applying the formulas of the previous subsection to the potential ( 3.20| ), a value of m ~ 10 13 GeV 
is required in order to obtain the observed amplitude of density fluctuations. Hence, the configura- 
tion space range with \(p\ > m p / but V((p) < m pl dominates the measure of field values. It can also 
be verified that the slow -roll trajectory is a local attractor in field initial condition space [|75|], even 
including metric fluctuations at the perturbative level [84]. 

With two scalar fields it is possible to construct a class of models which combine some of 
the nice features of large-field inflation (large phase space of initial conditions yielding inflation) 
and of small-field inflation (better contact with conventional particle physics). These are models 
of hybrid inflation [p5h. To give a prototypical example, consider two scalar fields (p and % with a 
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potential 

V(<p,z) = ^W-a 2 ) 2 + ^«V-ig 2 <py- (3-22) 

In the absence of thermal equilibrium, it is natural to assume that |<p| begins at large values, values 
for which the effective mass of % is positive and hence % begins at % = 0. The parameters in the 
potential ( 3.22| ) are now chosen such that <p is slowly rolling for values of |<p| somewhat smaller 



than m,pi, but that the potential energy for these field values is dominated by the first term on the 



right-hand side of ( |3.22[ ). The reader can easily verify that for this model it is no longer required 
to have values of |<p| greater than m p i in order to obtain slow-rolling The field (p is slowly rolling 
whereas the potential energy is determined by the contribution from X- Once |<p| drops to the value 

or 



(p c \ = ^-^o (3.23) 
8 

the configuration % = becomes unstable and decays to its ground state \% \ = a, yielding a graceful 
exit from inflation. Since in this example the ground state of % is not unique, there is the possibility 



of the formation of topological defects at the end of inflation (see [86, 87, 88] for reviews of 
topological defects in cosmology). 

3.4 Reheating in Inflationary Cosmology 

During the inflationary period the energy density of matter becomes completely dominated by 
the energy of the inflaton field. The density of regular matter (the kind of matter which is observed 
today) is exponentially diluted. An essential part of any inflationary model is the mechanism by 
which the energy density is transferred to regular matter at the end of the period of inflation. This is 
the so-called "reheating" process. It typically begins with a non-perturbative decay process called 
"preheating". We will give a brief overview of reheating in this subsection (for an in-depth recent 



review the reader is referred to [89] 



3.4.1 Preheating 

After the slow -roll conditions break down, the period of inflation ends, and the inflaton (p be- 
gins to oscillate around its ground state. Due to couplings of cp to other matter fields, the energy of 
the universe, which at the end of the period of inflation is stored completely in <p, gets transferred to 
the matter fields of the particle physics Standard Model. Initially, the energy transfer was described 
perturbatively [p0|]. Later, it was realized [ plj , p2| p3|, p4[ ] that through a parametric resonance in- 
stability, particles are very rapidly produced, leading to a fast energy transfer ("preheating"). The 
quanta later thermalize, and thereafter the universe evolves as described by SBB cosmology. 

Let us give a brief overview of the theory of reheating in inflationary cosmology. We assume 
that the inflaton (p is coupled to another scalar field % an d take the interaction Lagrangian to be 

^nt = -^V* 2 , (3.24) 
where g is a dimensionless coupling constant. 



8 Note that the slow-roll conditions ( p.lq ) and ( 3.17 ) were derived assuming that H is given by the contribution of 
<p to V which is not the case here. 
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In the initial perturbative analysis of reheating [[Xj the decay of the coherently oscillating 
inflaton field cp is treated perturbatively. The interaction rate T of processes in which two cp quanta 
interact to produce a pair of % particles is taken as the decay rate of the inflaton. The inflaton 
dynamics is modeled via an effective equation 

(p + 3H(p + r<p = -V\(p). (3.25) 

For small coupling constant, the interaction rate T is typically much smaller than the Hubble pa- 
rameter at the end of inflation. Thus, at the beginning of the phase of inflaton oscillations, the 
energy loss into particles is initially negligible compared to the energy loss due to the expansion 
of space. It is only once the Hubble expansion rate decreases to a value comparable to T that % 
particle production becomes effective. It is the energy density at the time when H = Y which de- 
termines how much energy ends up in % particles and thus determines the "reheating temperature", 
the temperature of the SM fields after energy transfer. 

T R ~ {rm p i) l/2 . (3.26) 

Since T is proportional to the square of the coupling constant g which is generally very small, per- 
turbative reheating is slow and produces a reheating temperature which can be very low compared 
to the energy scale at which inflation takes place. 

There are two main problems with the perturbative decay analysis described above. First of 
all, even if the inflaton decay were perturbative, it is not justified to use the heuristic equation ( |3.25 ) 



since it violates the fluctuation-dissipation theorem: in systems with dissipation, there are always 



fluctuations, and these are missing in (3.25). For an improved effective equation of motion see e.g. 



The main problem with the perturbative analysis is that it does not take into account the co- 
herent nature of the inflaton field (p. At the beginning of the period of oscillations cp is not a 
superposition of free asymptotic single inflaton states, but rather a coherently oscillating homoge- 
neous field. The large amplitude of oscillation implies that it is well justified to treat the inflaton 
classically. However, the matter fields can be assumed to start off in their vacuum state (the red- 
shifting during the period of inflation will remove any matter particles present at the beginning of 
inflation). Thus, matter fields % must be treated quantum mechanically. The improved approach 



to reheating initiated in Q9 1[ | (see also [Pq]) is to consider reheating as a quantum production of % 
particles in a classical (p background. 

We expand the field % is terms of the usual creation and annihilation operators. The mode 
functions satisfy the equation 

Zk + 3Hz+(^+m 2 x + g 2 ^(t) 2 sin 2 {mt^jXk = 0, (3.27) 

where <t>(t) is the amplitude of oscillation of (p. 

The first level of approximation [^l| is to neglect the expansion of space, to estimate the 
efficiency of reheating in this approximation, and to check self-consistency of the approximation. 



The equation (3.27) then reduces to 



Xk + ( k 2 + m\ + gW sin 2 (mt))xk = 0, (3.28) 
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and <!>(?) is constant. The equation is the Mathieu equation [97] which is conventionally written in 
the form 

Xk + (A* - 2<?cos 2z)Xk = , (3.29) 

where we have introduced the dimensionless time variable z = mt and a prime now denotes the 
derivative with respect to z. Comparing the coefficients, we see that 

k 2 + m\ o 2 d> 2 

A k = T A + 2q q = ~r^- (3-30) 

ra z 4m z 



The equation ( p.29| ) describes parametric resonance in classical mechanics. From the examples 
in classical mechanics in which this equation arises we know that there are bands of k values 
(depending on the frequency m of the external source) which exhibit exponential instability. For 
the general theory of the Mathieu equation the reader is referred to [p7|]. For values of k in a 
resonance band, the growth of Xk can be written as 

Xk~e^, (3.31) 

where \i k is called the Floquet exponent. In the model we are considering in this subsection, 
resonance occurs for all long wavelength modes - roughly speaking all modes with 

k 2 <.gm<S>. (3.32) 

This is called broad-band resonance [p2,|94]. The Floquet exponent is of order 1, and this implies 
very efficient energy transfer from the coherently oscillating inflaton field to a gas of x particles. 
This initial stage of energy transfer is called preheating. The time scale of the energy transfer 



is short compared to the cosmological time scale H at the end of inflation [91]. Hence, the 
approximation of neglecting the expansion of space is self-consistent. 

It is possible to perform an improved analysis which keeps track of the expansion of space. 
The first step [ p4| ] is to rescale the field variable to extract the cosmological red-shifting 

X k (t)=a(t)x k (t), (3.33) 
and to work in terms of conformal time. The field X k obeys an equation without cosmological 



damping term which is very similar to ( p.29| ), except that the oscillatory correction term to the 
mass is replaced by a more general periodic function of the conformal time v\ . The Floquet theory 
of [ f)7| ] also applies to the resulting equation, and we reach the conclusion that the exponential 
instability of x (modulated by a(t)) persists when including the effects of the expansion of space. 
For details the reader is referred to [p4|] or to the recent review [J$S|]. 



3.4.2 Preheating of Metric Perturbations 

A parametric resonance instability should be expected for all fields which couple to the oscil- 
lating inflaton condensate. In particular, metric fluctuations should also be effected. As was shown 



in [54], there is no parametric resonance instability for long wavelength modes in the case of purely 



adiabatic perturbations (see also [38]). However, in the presence of entropy fluctuations, parametric 



amplification of cosmological perturbations driven by the oscillating inflaton condensate is possible 



7, 58] (see also [B8D) 
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To study the preheating of metric fluctuations, we refer back to the formalism of cosmological 



perturbations developed in Subsection ( |2.3.4| ). In the presence of an entropy perturbation 8S, the 



"conservation" equation ( [2.5 1[ ) for the variable £ (which describes the curvature perturbation in 
comoving gauge) gets replaced by 

P 



C 



P + P 



-8S. 



(3.34) 



If long wavelength fluctuations of a matter field % which coupled to the inflaton q> (% corresponds 
to an entropy mode) undergo exponential growth during reheating as they do in the case of a model 
with broad parametric resonance, then this will induce an exponential growth of 8S which will via 



(3.34) induce an exponentially growing curvature perturbation. 



For an elegant formalism to compute the source term in ( |3.34| ) in terms of the scalar fields q> 



and x the reader is referred to [p9|. This resonant growth of entropy fluctuations is only important 
in models in which the entropy fluctuations are not suppressed during inflation. Some recent ex- 



amples were studied in Q100|, [lOjJ , |102j|. The source for this instability need not be the oscillations 
of the inflaton field. In SUSY models, the decay of flat directions can also induce this instability 



[103]. 



3.5 Problems of Inflation 

In spite of the phenomenological success of the inflationary paradigm, conventional scalar 
field-driven inflation suffers from several important conceptual problems. 

The first problem concern the nature of the inflaton, the scalar field which generates the in- 
flationary expansion. No particle corresponding to the excitation of a scalar field has yet been 
observed in nature, and the Higgs field which is introduced to give elementary particles masses in 
the Standard Model of particle physics does not have the required flatness of the potential to yield 



inflation, unless it is non-minimally coupled to gravity [104]. In particle physics theories beyond 
the Standard Model there are often many scalar fields, but it is in general very hard to obtain the 
required flatness properties on the potential 

The second problem (the amplitude problem) relates to the amplitude of the spectrum of 
cosmological perturbations. In a wide class of inflationary models, obtaining the correct amplitude 



requires the introduction of a hierarchy in scales, namely [105] 



V(9) 
A<p 4 



< 10 



-12 



(3.35) 



where A<p is the change in the inflaton field during the minimal length of the inflationary period, 
and V((p) is the potential energy during inflation. 



A more serious problem is the trans-Planckian problem [ jlQq ]. Returning to the space-time 



diagram of inflation (see Figure |10J), we can immediately deduce that, provided that the period 
of inflation lasted sufficiently long (for GUT scale inflation the number is about 70 e-foldings), 
then all scales inside the Hubble radius today started out with a physical wavelength smaller than 
the Planck scale at the beginning of inflation. Now, the theory of cosmological perturbations is 
based on Einstein's theory of General Relativity coupled to a simple semi-classical description of 
matter. It is clear that these building blocks of the theory are inapplicable on scales comparable and 
smaller than the Planck scale. Thus, the key successful prediction of inflation (the theory of the 
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Figure 10: Space-time diagram (sketch) of inflationary cosmology where we have added an extra length 
scale, namely the Planck length l p i (majenta vertical line). The symbols have the same meaning as in 
Figure 5. Note, specifically, that - as long as the period of inflation lasts a couple of e-foldings longer 
than the minimal value required for inflation to address the problems of standard big bang cosmology - all 
wavelengths of cosmological interest to us today start out at the beginning of the period of inflation with a 
wavelength which is smaller than the Planck length. 



origin of fluctuations) is based on suspect calculations since new physics must enter into a correct 
computation of the spectrum of cosmological perturbations. The key question is as to whether the 
predictions obtained using the current theory are sensitive to the specifics of the unknown theory 
which takes over on small scales. 

One approach to study the sensitivity of the usual predictions of inflationary cosmology to 
the unknown physics on trans-Planckian scales is to study toy models of ultraviolet physics which 
allow explicit calculations. The first approach which was used [107, 108] is to replace the usual lin- 
ear dispersion relation for the Fourier modes of the fluctuations by a modified dispersion relation, a 
dispersion relation which is linear for physical wavenumbers smaller than the scale of new physics, 
but deviates on larger scales. Such dispersion relations were used previously to test the sensitivity 
of black hole radiation on the unknown physics of the UV [109, 110]. It was found [107] that if 
the evolution of modes on the trans-Planckian scales is non-adiabatic, then substantial deviations 
of the spectrum of fluctuations from the usual results are possible. Non-adiabatic evolution turns 
an initial state minimizing the energy density into a state which is excited once the wavelength 
becomes larger than the cutoff scale. Back-reaction effects of these excitations may limit the mag- 
nitude of the trans-Planckian effects if we assume that inflation took place [ 111 , 112 , 1 13 ]. On the 
other hand, large trans-Planckian effects may also prevent the onset of inflation. 

A fourth problem is the singularity problem. It was known for a long time that standard Big 
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Bang cosmology cannot be the complete story of the early universe because of the initial singularity, 
a singularity which is unavoidable when basing cosmology on Einstein's field equations in the 



presence of a matter source obeying the weak energy conditions (see e.g. [ |114j ] for a textbook 
discussion). Recently, the singularity theorems have been generalized to apply to Einstein gravity 
coupled to scalar field matter, i.e. to scalar field-driven inflationary cosmology [ ]115| ]. It was shown 
that, in this context, a past singularity at some point in space is unavoidable. Thus we know, from 
the outset, that scalar field-driven inflation cannot be the ultimate theory of the very early universe. 

The Achilles heel of scalar field-driven inflationary cosmology may be the cosmological con- 
stant problem. We know from observations that the large quantum vacuum energy of field theo- 
ries does not gravitate today. However, to obtain a period of inflation one is using the part of the 
energy-momentum tensor of the scalar field which looks like the vacuum energy. In the absence of 
a solution of the cosmological constant problem it is unclear whether scalar field-driven inflation is 
robust, i.e. whether the mechanism which renders the quantum vacuum energy gravitationally inert 
today will not also prevent the vacuum energy from gravitating during the period of slow-rolling of 
the inflaton field. Note that the approach to addressing the cosmological constant problem making 



use of the gravitational back-reaction of long range fluctuations (see [ 116] for a summary of this 
approach) does not prevent a long period of inflation in the early universe. 

A final problem which we will mention here is the concern that the energy scale at which 
inflation takes place is too high to justify an effective field theory analysis based on Einstein gravity. 
In simple toy models of inflation, the energy scale during the period of inflation is about 10 16 GeV, 
very close to the string scale in many string models, and not too far from the Planck scale. Thus, 
correction terms in the effective action for matter and gravity may already be important at the 
energy scale of inflation, and the cosmological dynamics may be rather different from what is 
obtained when neglecting the correction terms. 

In Figure pi] we show once again the space-time sketch of inflationary cosmology. In addition 
to the length scales which appear in the previous versions of this figure, we have now shaded the 
"zones of ignorance", zones where the Einstein gravity effective action is sure to break down. As 
described above, fluctuations emerge from the short distance zone of ignorance (except if the period 
of inflation is very short), and the energy scale of inflation might put the period of inflation too close 
to the high energy density zone of ignorance to trust the predictions based on using the Einstein 
action. 

The arguments in this subsection provide a motivation for considering alternative scenarios of 
early universe cosmology. Below we will focus on two scenarios, the matter bounce and string gas 
cosmology. It is important to emphasize that these are not the only alternative scenarios. Others 



include the Pre-Big-Bang scenario [117] and the Ekpyrotic paradigm [118]. 



4. Matter Bounce 

The matter bounce scenario is a non-singular cosmology in which time runs from — oo to +oo. 
Negative times correspond to a contracting phase, positive times to expansion. It is assumed that 
the bounce phase is short 9 . 



Short means short compared to the time it takes light to travel over distances of about 1mm, the physical length 
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Figure 11: Space-time diagram (sketch) of inflationary cosmology including the two zones of ignorance - 
sub-Planckian wavelengths and trans-Planckian densities. The symbols have the same meaning as in Figure 
5. Note, specifically, that - as long as the period of inflation lasts a couple of e-foldings longer than the min- 
imal value required for inflation to address the problems of standard big bang cosmology - all wavelengths 
of cosmological interest to us today start out at the beginning of the period of inflation with a wavelength 
which is in the zone of ignorance. 



As indicated in Figure ^, the fluctuations which we observe today have exited the Hubble 
radius early in the contracting phase. The word matter in matter bounce should emphasize that we 
are assuming that the phase during which the fluctuations which we observe today are exiting the 
Hubble radius is dominated by cold pressure-less matter. Viewing the contracting phase as the time 
reverse of the expanding phase we are living in, this last assumption we have made does not seem 
very restrictive. In fact, one could expect that some entropy is generated during the bounce, in 
which case the matter phase during the contracting period would last up to higher energy densities 
than it does in the expanding phase. 

4.1 Models for a Matter Bounce and Background Cosmology 



It follows from the singularity theorems of General Relativity (see e.g. [|114J]) that in order 
to obtain a non-singular bounce we must either invoke matter which violates the weak energy 
condition, or go beyond General Relativity. A review of ways to obtain a non-singular cosmology 



is given in [20]. The major challenges in constructing a non-singular bouncing cosmology include 



which corresponds to to the largest comoving scales we observe today (assuming that the energy scale of the bounce is 
about 10 16 GeV). 
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the following: first, the model should be free of ghost-like excitations, since a model with ghosts 
will be unstable [119]. Next, the background cosmology should be stable towards the addition of 
radiation to the matter sector. Finally, the homogeneous and isotropic background trajectory should 
also be stable against anisotropic stress. 

In the following we will only mention a few recent attempts which the author of these lecture 
notes has been involved in. We will first discuss models obtained by modifying the matter sector 
of the theory, and then models where the gravitational sector of the theory is changed. 

Introducing quintom matter [ ]120| ] yields a way of obtaining a non-singular bouncing cosmol- 
ogy, as discussed in [p^]. Quintom matter is a set of two matter fields, one of them regular matter 
(obeying the weak energy condition), the second a field with opposite sign kinetic term, a field 
which violates the energy conditions. We can fll9[ ] model both matter components with scalar 
fields. Let us denote the mass of the regular one ((p) by m, and by M that of the field (p with wrong 
sign kinetic term. We assume that early in the contracting phase both fields are oscillating, but that 
the amplitude stf of (p greatly exceeds the corresponding amplitude of (p such that the energy 
density is dominated by (p. Both fields will initially be oscillating during the contracting phase, 
and both amplitudes grow at the same rate. At some point, srf will become so large that the oscil- 
lations of cp freeze out Then, 80 will grow only slowly, whereas s$ will continue to increase. 
Thus, the (negative) energy density in <p will grow in absolute values relative to that of (p. The total 
energy density will decrease towards 0. At that point, H = by the Friedmann equations. It can 
in fact easily be seen that H > when H = 0. Hence, a non-singular bounce occurs. The Higgs 
sector of the Lee-Wick model [121 ] provides a concrete realization of the quintom bounce model, 
as studied in [|25|]. Quintom models like all other models with negative sign kinetic terms suffer 



from an instability problem [ |119| ] in the matter sector and are hence problematic. In addition, they 
are unstable against the addition of radiation (see e.g. [122]) and anisotropic stress. 

An improved way of obtaining a non-singular bouncing cosmology using modified matter is 
by using a ghost condensate field [123]. The ghost condensation mechanism is the analog of the 
Higgs mechanism in the kinetic sector of the theory. In the Higgs mechanism we take a field (j) 
whose mass when evaluated at = is tachyonic, add higher powers of (p 2 to the potential term 
in the Lagrangian such that there is a stable fixed point = v ^ 0, and thus when expanded about 
= v the mass term has the "safe" non-tachyonic sign. In the ghost condensate construction we 
take a field whose kinetic term 

X = -g^drfd^ (4.1) 

appears with the wrong sign in the Lagrangian. Then, we add higher powers of X to the kinetic 
Lagrangian such that there is a stable fixed point X = c 2 and such that when expanded about X = c 2 
the fluctuations have the regular sign of the kinetic term: 

S£= l -M\X-c 2 f-V (4.2) 

where V (<p) is a usual potential function, M is a characteristic mass scale and the dimensions of 
are chosen such that X is dimensionless. 

In the context of cosmology, the ghost condensate is 

= ct (4.3) 



10' 



This corresponds to the time reverse of entering a region of large-field inflation. 
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and breaks local Lorentz invariance. Now let us expand the homogeneous component of about 
the ghost condensate: 

<j)(t) =ct + n(t). (4.4) 

If ft < then the gravitational energy density is negative, and a non-singular bounce is possible. 
Thus, in [ |123| ] we constructed a model in which the ghost condensate field starts at negative values 
and the potential V(<p) is negligible. As <p approaches <p = it encounters a positive potential 
which slows it down, leading to iz < and hence to negative gravitational energy density. Thus, a 
non-singular bounce can occur. We take the potential to be of the form 

V(0) ~ <T a (4-5) 

for 1 1 S> M, where M is the mass scale above which the higher derivative kinetic terms are impor- 
tant. For sufficiently large values of a, namely 

a > 6 , (4.6) 

the energy density in the ghost condensate increases faster than that of radiation and anisotropic 
stress at the universe contracts . Hence, this bouncing cosmology is stable against the addition of 
radiation and anisotropic stress. 

Turning to the second way of obtaining a non-singular bounce, namely by modifying the grav- 
itational action, we should emphasize that modifications of the gravitational action are expected at 
the high densities at which the bounce will occur. Concrete examples were already mentioned in the 
Introduction: the higher derivative Lagrangian resulting from the nonsingular universe construction 



of 017Q, the model of [16] which is based on a non-local higher derivative action which is ghost- 



free about Minkowski space-time, mirage cosmologies Q125| ] resulting from the effective action of 
gravity on a brane which is moving into and out of a high-curvature throat in a higher-dimensional 
space-time. 

More recently, interest has focused on a non-singular bouncing cosmology which emerges 



from Horava-Lifshitz gravity Q22[] . Horava-Lifshitz gravity is a power-counting renormalizable 
quantum theory of gravity which is based on anisotropic scaling between space and time. In one 
sense it is a very conservative approach to quantum gravity in the sense that the dynamical degrees 
of freedom are the usual ones: the spatial metric the lapse function N and the shift vector TV;. 
There are no extra dimensions or new degrees of freedom such as strings. On the other hand, it 
is a very radical approach in the sense that it gives up local Lorentz invariance and space-time 
diffeomorphism invariance. The residual symmetries are local rotational invariance and spatial 
diffeomorphism. On the other hand, there is a scaling symmetry under 

x ->■ Ix , t ->■ l 3 t , (4.7) 

where / is the scaling parameter. The Lagrangian is constructed to contain all terms which are 
consistent with the residual symmetries and with power-counting renormalizability with respect 
to the above scaling symmetry. The action does not contain any higher time derivative terms, 
but it does contain higher space derivative terms. These terms dominate in the ultraviolet. In the 
infrared, the extra terms are sub-dominant. Thus, in the infrared the theory is supposed to flow 
towards General Relativity. 
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The extra spatial derivative terms in the gravitational action have an important effect in the 
very early universe at high energy densities. As shown in [|2l|], if the spatial sections have a non- 
vanishing spatial curvature constant k, then the higher-derivative terms in the action will lead to 
terms in the Friedmann equations which act as ghost radiation and ghost anisotropic stress, i.e. 
terms of gravitational origin and negative effective energy density which scale as a~ 4 and a~ 6 , 
respectively. Starting with a contracting universe dominated by regular matter, eventually the ghost 
terms will catch up and yield a non-singular bounce in analogy to how the ghost matter in the 
quintom model does. 

The analysis of the spectrum of cosmological perturbations on scales of current observational 
interest is, however, independent of the details of the bouncing phase, as long as that phase is short 
compared to the time it takes light to travel over the length scales of current interest, and as long 
as no extra degrees of freedom (extra when compared to those arising in Einstein gravity) become 
important 1 1 . We now turn to a discussion of the evolution of fluctuations in a generic matter bounce 
model. 

4.2 Fluctuations in a Matter Bounce 

First we will consider fluctuations in a matter bounce without extra degrees of freedom. In this 
case, we need only focus on the usual fluctuation variable v. The equation of motion its Fourier 
mode Vk is 



, ( 2_Z_ 

Z 

If the equation of state of the background is time-independent, then z ~ a and hence the negative 



v k +{k 2 -j)v k = 0. (4.8) 



square mass term in ( |4.8[ ) is H . Thus, on length scales smaller than the Hubble radius, the solutions 
of (-k8) are oscillating, whereas on larger scales they are frozen in, and their amplitude depends on 
the time evolution of z. 

In the case of an expanding universe the dominant mode of v scales as z. However, in a 
contracting universe it is the second of the two modes which dominates. If the contracting phase is 
matter-dominated, i.e. a(t) ~ f 2 / 3 and r\(t) ~ f 1 / 3 the dominant mode of v scales as r\~ l and hence 

Vfc(ij) = ciTj 2 + c 2 ir\ (4-9) 

where c\ and ci are constants. The c\ mode is the mode for which £ is constant on super-Hubble 
scales. However, in a contracting universe it is the C2 mode which dominates and leads to a scale- 



invariant spectrum [ |13| , NJ , |1 



P f (*,Tj) ~ * 3 h(7?)|V 2 (T]) (4-10) 

-j^lv^C*))! 2 ^) 2 -* 3 - 1 - 2 



const . 



where we have made use of the scaling of the dominant mode of v k , the Hubble radius crossing 
condition ~ k~ l , and the assumption that we have a vacuum spectrum at Hubble radius 

crossing. 



1 1 This last condition is very relevant in the case of the Horava-Lifshitz bounce where we do have to worry about 
extra gravitational degrees of freedom. We will come back to this issue at the end of the following subsection. 
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Up to this point, the analysis shows that in the contracting phase the curvature fluctuations 
are scale-invariant. In non-singular bouncing cosmologies, the fluctuations can be followed in an 
unambiguous way through the bounce. This was done in the case of the higher derivative bounce 
model of fli6| ] in [23], and in a non-singular bouncing mirage cosmology model the analysis was 



performed in [18 ]. In the case of the Quintom and Lee-Wick bounces, respectively, the fluctuations 
were followed through the bounce in [24] and [25], respectively. Finally, in the case of the non- 
projectable version 12 of the HL bounce, the evolution of the fluctuations through the bounce was 
recently studied in [26]. 

The equations of motion can be solved numerically without approximation. Alternatively, 
we can solve them approximately using analytical techniques. Key to the analytical analysis are 
the General Relativistic matching conditions for fluctuations across a phase transition in the back- 
ground [126, 127]. These conditions imply that both <J> and £ are conserved at the bounce, where 



1 k 2 <S> 



(4.11) 



However, as stressed in [128], these matching conditions can only be used at a transition when the 
background metric obeys the matching conditions. This is not the case if we were to match directly 
between the contracting matter phase and the expanding matter phase. 

In the case of a non-singular bounce we have three phases: the initial contracting phase with 
a fixed equation of state (e.g. w = 0), a bounce phase during which the universe smoothly transits 
between contraction and expansion, and finally the expanding phase with constant w. We need 
to apply the matching conditions twice: first at the transition between the contracting phase and 
the bounce phase (on both sides of the matching surface the universe is contracting), and then 
between the bouncing phase and the expanding phase. The bottom line of the studies of [23, 18, 
24, 25, 26] is that on length scales large compared to the time of the bounce, the spectrum of 
curvature fluctuations is not changed during the bounce phase. Since typically the bounce time 
is set by a microphysical scale whereas the wavelength of fluctuations which we observe today is 
macroscopic (about 1mm if the bounce scale is set by the particle physics GUT scale), we conclude 
that for scales relevant to current observations the spectrum is unchanged during the bounce. This 
completes the demonstration that a non-singular matter bounce leads to a scale-invariant spectrum 
of cosmological perturbations after the bounce provided that the initial spectrum on sub-Hubble 
scales is vacuum. Initial thermal fluctuations were followed through the bounce in [27]. Note that 
perturbations are processed during a bounce and this implies that even if the background cosmology 
were cyclic, the full evolution including linear fluctuations would be non-cyclic [32]. 

The above analysis is applicable only as long as no new degrees of freedom become relevant 
at high energy densities, in particular during the bounce phase. In non-singular bounce models 
obtained by modifying the matter sector, new degrees of freedom arise from the extra matter fields. 
They can thus give entropy fluctuations which may compete with the adiabatic mode studied above. 
In the quintom bounce model this issue has recently been studied in [129]. It was found that 
fluctuations in the ghost field which yields the bounce are unimportant on large scales since they 
have a blue spectrum. However, entropy fluctuations due to extra low-mass fields can be important. 
They yield the "matter bounce curvaton" mechanism. Their spectrum is also scale-invariant. 



12 



See below for the specification of this case. 
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In non-singular bouncing models obtained by modifying the gravitational sector of the theory 
the identification of potential extra degrees of freedom is more difficult. As an example, let us 
mention the situation in the case of the Horava-Lifshitz bounce. The theory has the same number 
of geometric degrees of freedom as General Relativity, but less symmetries. Thus, more of the 
degrees of freedom are physical. Recall from the discussion of the theory of cosmological pertur- 
bations in Section 2 that there are ten total geometrical degrees of freedom for linear cosmological 
perturbations, four of them being scalar, four vector and two tensor. In Einstein gravity the sym- 
metry group of space-time diffeomorphisms is generated at the level of linear fluctuations by four 
functions, leaving six of the ten geometrical variables as physical - two scalar, two vector and two 
tensor modes. In the absence of anisotropic stress the number of scalar variables is reduced by one, 
and the Hamiltonian constraint relates the remaining scalar metric fluctuation to matter. 

In Horava-Lifshitz gravity one loses one scalar gauge degree of freedom, namely that of space- 
dependent time reparametrizations. Thus, one expects an extra physical degree of freedom. It has 



been recently been shown [130] that in the projectable version of the theory (in which the lapse 
function N(t) is constrainted to be a function of time only) the extra degree of scalar cosmolog- 
ical perturbations is either ghost-like or tachyonic, depending on parameters in the Lagrangian. 
Thus, the theory appears to be ill-behaved in the context of cosmology. However, in the full non- 
projectable version (in which the lapse N(t,x) is a function of both space and time, the extra degree 



of freedom is well behaved. It is important on ultraviolet scales but decouples in the infrared [131] 



4.3 Key Predictions of the Matter Bounce 

Canonical single field inflation models predict very small non-Gaussianities in the spectrum 
of fluctuations. One way to characterize the non-Gaussianities is via the three point function of the 



curvature fluctuation, also called the "bispectrum". As realized in [ |132| ], the bispectrum induced 
in the minimal matter bounce scenario (no entropy modes considered) has an amplitude which is 
at the borderline of what the Planck satellite experiment will be able to detect, and it has a special 
form. These are specific predictions of the matter bounce scenario with which the matter bounce 



scenario can be distinguished from those of standard inflationary models (see [133] for a recent 
detailed review of non-Gaussianities in inflationary cosmology and a list of references). In the 
following we give a very brief summary of the analysis of non-Gaussianities in the matter bounce 
scenario. 

Non-Gaussianities are induced in any cosmological model simply because the Einstein equa- 
tions are non-linear. In momentum space, the bispectrum contains amplitude and shape informa- 
tion. The bispectrum is a function of the three momenta. Momentum conservation implies that the 
three momenta have to add up to zero. However, this still leaves a rich shape information in the 
bispectrum in addition to the information about the overall amplitude. 

A formalism to compute the non-Gaussianities for the curvature fluctuation variable £ was 



developed in [ ]134| ]. Working in the interaction representation, the three-point function of £ is 



given to leading order by 

< C(f,*i)C(f,fc)C(f,* 3 ) > (4.12) 

= if'dt' < [C(/,*i)Cfrfc)CM 3 ),i**(0] > » 
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where ? ; corresponds to the initial time before which there are any non-Gaussianities. The square 
parentheses indicate the commutator, and L[ nt is the interaction Lagrangian 

The interaction Lagrangian contains many terms. In particular, there are terms containing the 
time derivative of £. Each term leads to a particular shape of the bispectrum. In an expanding uni- 
verse such as in inflationary cosmology £ = 0. However, in a contracting phase the time derivative 
of £ does not vanish since the dominant mode is growing in time. Hence, there are new contribu- 
tions to the shape which have a very different form from the shape of the terms which appear in 
inflationary cosmology. The larger value of the amplitude of the bispectrum follows again from the 
fact that there is a mode function which grows in time in the contracting phase. 

The three-point function can be expressed in the following general form: 

p2 

< C(*i)C(fe)C(*3) > = (2x) 7 8(£l" " 



, (4.13) 

where ki = and si is the shape function. In this expression we have factored out the dependence 
on the power spectrum Pr. In inflationary cosmology it has become usual to express the bispectrum 
in terms of a non-Gaussianity parameter fyi- However, this is only useful if the shape of the three 
point function is known. As a generalization, we here use [|132|] 



m.SuhM) = g ^W) . (4.14, 

The computation of the bispectrum is tedious. In the case of the matter bounce (no entropy 
fluctuations) the result is 



256rm u k J 



¥J 

-9£*f*5 + 5£*fjfcj (4.15) 

¥J ¥J 



. £ kfkjkl + 9 £ kfkjk, 



This equation describes the shape which is predicted. Some of the terms (such as the last two) are 
the same as those which occur in single field slow-roll inflation, but the others are new. Note, in 
particular, that the new terms are not negligible. 

If we project the resulting shape function srf onto some popular shape masks we 



tt = - 3 4 > ( 4 - 16 ) 



for the local shape (&i <C &2 = k^). This is negative and of order 0(1). For the equilateral form 
(&i =k2 = ki) the result is 



I equil 



255 
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and for the folded form {k\ = 2&2 = 2^3) one obtains the value 



IT" = -\ • (4-18) 



These amplitudes are close to what the Planck CMB satellite experiment will be able to detect. 
The matter bounce scenario also predicts a change in the slope of the primordial power spec- 



trum on small scales [ |135[ ] : scales which exit the Hubble radius in the radiation phase retain a blue 
spectrum since the squeezing rate on scales larger than the Hubble radius is insufficient to give 
longer wavelength modes a sufficient boost relative to the shorter wavelength ones. 

4.4 Problems of the Matter Bounce 

The main problem of the matter bounce scenario is that the physics which yields the non- 
singular bounce is not (yet) well established. It is clearly required to go beyond conventional 
quantum field theory coupled to General Relativity to obtain such a bounce. New Planck-scale 
physics needs to be invoked to obtain the required background cosmology. We remind the reader 
that it is precisely the absence of such new physics which needs to be invoked to argue for an 
inflationary background evolution. 

In terms of the trans-Planckian problem for cosmological fluctuations, the bouncing cosmol- 
ogy has a clear advantage: the wavelength of the fluctuations which we are interested in always 
remains in the far infrared compared to the ultraviolet scales of the new physics. It can be shown 
that the correction terms due to the new ultraviolet physics on the evolution of fluctuations of inter- 



est to observers are negligible (see e.g. [26]). 

The second main problem of the matter bounce scenario may be the sensitivity of the bounce 
to assumptions on the initial conditions in the far past. The amplitude of the classical initial fluctu- 
ations must be small in order for the homogeneous background dynamics not to be perturbed, and 



for the vacuum fluctuations to dominate on the scales relevant in the infrared (see e.g. [136] for 
a criticism of initial conditions in a related scenario, the Pre-Big-Bang scenario). In addition, the 
initial shear must be very small. In some bouncing cosmology backgrounds the bouncing solution 
is unstable to the smallest addition of shear. This, however, is not the case in other models such as 



the Horava-Lifshitz bounce [ |21[ ] or the ghost condensate bounce [ |123| ]. 

The matter bounce scenario does not address the cosmological constant problem, but on the 
other hand the scenario does not depend on how the cosmological constant problem is solved. Note 
that inflationary cosmology is NOT robust in this respect. 



5. String Gas Cosmology 

Another scenario of early universe cosmology which can provide an explanation for the current 



data is string gas cosmology. It is an "emergent universe" [137] scenario in which the universe 
begins in a long hot and almost static phase. The key input from string theory is that matter is a gas 
of closed fundamental strings compared to a gas of point particles as is assumed in Standard Big 
Bang cosmology. 
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5.1 Principles and Background 
5.1.1 Principles 

String theory may be the best candidate we have at the present time for a quantum theory 
of gravity which unifies all four forces of nature. This theory is, however, currently not yet well 
understood beyond perturbation theory 13 . For applications to early universe cosmology, however, 
a non-perturbative understanding will be essential. 

In the absence of a non-perturbative formulation of string theory, the approach to string cos- 
mology which we have suggested, string gas cosmology [ |2g| ] (see also p9|], and [11], |139| , |140fl for 



reviews and more complete list of references), is to focus on symmetries and degrees of freedom 
which are new to string theory (compared to point particle theories) and which will be part of any 
non-perturbative string theory, and to use them to develop a new cosmology. The symmetry we 
make use of is T-duality, and the new degrees of freedom are the string oscillatory modes and the 
string winding modes. 

String gas cosmology is based on coupling a classical background which includes the graviton 
and the dilaton fields to a gas of closed strings (and possibly other basic degrees of freedom of 



string theory such as "branes" [ ]141p ). All dimensions of space are taken to be compact, for reasons 
which will become clear later. For simplicity, we take all spatial directions to be toroidal and denote 
the radius of the torus by R. Strings have three types of states: momentum modes which represent 
the center of mass motion of the string, oscillatory modes which represent the fluctuations of the 
strings, and winding modes counting the number of times a string wraps the torus. 

Since the number of string oscillatory states increases exponentially with energy, there is a 



limiting temperature for a gas of strings in thermal equilibrium, the Hagedorn temperature [3C] 
Th- Thus, if we take a box of strings and adiabatically decrease the box size, the temperature 
will never diverge. This is the first indication that string theory has the potential to resolve the 
cosmological singularity problem. 

The second key feature of string theory upon which string gas cosmology is based is T-duality. 
To introduce this symmetry, let us discuss the radius dependence of the energy of the basic string 
states: The energy of an oscillatory mode is independent of R, momentum mode energies are 
quantized in units of l/R, i.e. 

E n = n^, (5.1) 
K 

where l s is the suing length and pt is the mass per unit length of a string. The winding mode 
energies are quantized in units of R, i.e. 

E m = m/j,R, (5.2) 

where both n and m are integers. Thus, a new symmetry of the spectrum of string states emerges: 
Under the change 

R -> l/R (5.3) 



Note that there are concrete proposals for a non-perturbative formulation such as the AdS/CFT correspondence 

pi. 
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in the radius of the torus (in units of l s ) the energy spectrum of string states is invariant if winding 
and momentum quantum numbers are interchanged 

(n,m) — > (m,n) . (5.4) 

The above symmetry is the simplest element of a larger symmetry group, the T-duality symmetry 
group which in general also mixes fluxes and geometry. The string vertex operators are consistent 
with this symmetry, and thus T-duality is a symmetry of perturbative string theory. Postulating that 
T-duality extends to non-perturbative string theory leads [ ]142| ] to the need of adding D-branes to 
the list of fundamental objects in string theory. With this addition, T-duality is expected to be a 
symmetry of non-perturbative string theory. Specifically, T-duality will take a spectrum of stable 
Type IIA branes and map it into a corresponding spectrum of stable Type IIB branes with identical 



masses [143]. 



As discussed in \UM, the above T-duality symmetry leads to an equivalence between small and 



large spaces, an equivalence elaborated on further in [144, 145]. 



5.1.2 Background Cosmology 

That string gas cosmology will lead to a dynamical evolution of the early universe very differ- 
ent from what is obtained in standard and inflationary cosmology can already be seen by combining 
the basic ingredients from string theory discussed in the previous subsection. As the radius of a 
box of strings decreases from an initially very large value - maintaining thermal equilibrium - , the 
temperature first rises as in standard cosmology since the string states which are occupied (the mo- 
mentum modes) get heavier. However, as the temperature approaches the Hagedorn temperature, 
the energy begins to flow into the oscillatory modes and the increase in temperature levels off. As 
the radius R decreases below the string scale, the temperature begins to decrease as the energy be- 



gins to flow into the winding modes whose energy decreases as R decreases (see Figure |12J). Thus, 



as argued in H28Q, the temperature singularity of early universe cosmology should be resolved in 
string gas cosmology. 

The equations that govern the background of string gas cosmology are not known. The Ein- 
stein equations are not the correct equations since they do not obey the T-duality symmetry of string 
theory. Many early studies of string gas cosmology were based on using the dilaton gravity equa- 



tions [ ]146| , [147], |148| ]. However, these equations are not satisfactory, either. Firstly, we expect that 
string theory correction terms to the low energy effective action of string theory become dominant 
in the Hagedorn phase. Secondly, the dilaton gravity equations yield a rapidly changing dilaton 
during the Hagedorn phase (in the string frame). Once the dilaton becomes large, it becomes in- 
consistent to focus on fundamental string states rather than brane states. In other words, using 
dilaton gravity as a background for string gas cosmology does not correctly reflect the S-duality 
symmetry of string theory. Recently, a background for string gas cosmology including a rolling 



tachyon was proposed [ |149[ ] which allows a background in the Hagedorn phase with constant scale 
factor and constant dilaton. Another study of this problem was given in [ |150| ]. 

Some conclusions about the time-temperature relation in string gas cosmology can be derived 
based on thermodynamical considerations alone. One possibility is that R starts out much smaller 
than the self-dual value and increases monotonically. From Figure [12] it then follows that the time- 
temperature curve will correspond to that of a bouncing cosmology. Alternatively, it is possible 
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Figure 12: The temperature (vertical axis) as a function of radius (horizontal axis) of a gas of closed strings 
in thermal equilibrium. Note the absence of a temperature singularity. The range of values of R for which 
the temperature is close to the Hagedorn temperature Th depends on the total entropy of the universe. The 
upper of the two curves corresponds to a universe with larger entropy. 




> t 



p = R P = rho / 3 



Figure 13: The dynamics of string gas cosmology. The vertical axis represents the scale factor of the 
universe, the horizontal axis is time. Along the horizontal axis, the approximate equation of state is also 
indicated. During the Hagedorn phase the pressure is negligible due to the cancellation between the positive 
pressure of the momentum modes and the negative pressure of the winding modes, after time tR the equation 
of state is that of a radiation-dominated universe. 



that the universe starts out in a meta-stable state near the Hagedorn temperature, the Hagedorn 
phase, and then smoothly evolves into an expanding phase dominated by radiation like in standard 
cosmology (Figure 13). Note that we are assuming that not only the scale factor but also the dilaton 
is constant in time. 

The transition between the quasi-static Hagedorn phase and the radiation phase at the time f# is 
a consequence of the annihilation of string winding modes into string loops (see Figure 14). Since 
this process corresponds to the production of radiation, we denote this time by the same symbol 
tR as the time of reheating in inflationary cosmology. As pointed out in p8|], this annihilation 
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Figure 14: The process by which string loops are produced via the intersection of winding strings. The 
top and bottom lines are identified and the space between these lines represents space with one toroidal 
dimension un-wrapped. 

process only is possible in at most three large spatial dimensions. This is a simple dimension 
counting argument: string world sheets have measure zero intersection probability in more than 
four large space-time dimensions. Hence, string gas cosmology may provide a natural mechanism 
for explaining why there are exactly three large spatial dimensions. This argument was supported 



by numerical studies of string evolution in three and four spatial dimensions Q151| ] (see also [ |152[ ]). 



The flow of energy from winding modes to string loops can be modelled by effective Boltzmann 



equations [153] analogous to those used to describe the flow of energy between infinite cosmic 



strings and cosmic string loops (see e.g. [86, K7L K8h for reviews). 



Several comments are in place concerning the above mechanism. First, in the analysis of 
[153] it was assumed that the string interaction rates were time-independent. If the dynamics of 



the Hagedorn phase is modelled by dilaton gravity, the dilaton is rapidly changing during the phase 



in which the string frame scale factor is static. As discussed in [154, 155] (see also fl!56[ ]), in 
this case the mechanism which tells us that exactly three spatial dimensions become macroscopic 
does not work. Another comment concerns the isotropy of the three large dimensions. In contrast 
to the situation in Standard cosmology, in string gas cosmology the anisotropy decreases in the 



expanding phase [ ]157| ]. Thus, there is a natural isotropization mechanism for the three large spatial 
dimensions. 

5.1.3 Moduli Stabilization 

In the following, we shall assume that either the mechanism of ||2^] for setting in motion 
the preferential expansion of exactly three spatial dimensions works, or, alternatively, that three 
dimensions are distinguished from the beginning as being large. In either case, we must address 
the moduli stabilization problem, i.e. we must show that the radii (radions) and shapes of the extra 
dimensions are stabilized. This is a major challenge in string-motivated field theory models of 
higher dimensions. The situation in string gas cosmology is much better in comparison: all moduli 
except for the dilaton are stabilized without the need of introducing extra ingredients such as fluxes 
or special non-perturbative effects. 

Radion stabilization in the string frame was initially studied in [ ]158| ]. The basic idea is that 
the winding modes about the extra spatial dimensions provide a confining force which prevents the 
radii from increasing whereas the momentum modes provide a force which resists the complete 
contraction. Thus, there will be a stable minimum of the effective potential for the radion. 

In order to make contact with late time cosmology, it is important to consider the issue of 
radion stabilization when the dilaton is frozen, or, more generally, in the Einstein frame. As was 
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discussed in [159, 160] (see also earlier comments in [158]), the existence of string modes which 



are massless at the self-dual radius is crucial in obtaining radion stabilization in the Einstein frame 



(for more general studies of the importance of massless modes in string cosmology see Q161| , |162| ]). 



Such massless modes do not exist in all known string theories. They exist in the Heterotic theory, 



but not in Type II theories [ |142[ ]. The following discussion is taken from [ |163[ ]. 

Let us consider the equations of motion which arise from coupling the Einstein action to a 
string gas. In the anisotropic setting when the metric is taken to be 

6 

ds 1 = dt 2 - a(t) 2 dx 2 - £ b a {t) 2 dy 2 a , (5.5) 
a=l 

where the y a are the internal coordinates, the equation of motion for b a becomes 

b a + (3H+ £ b ±)b a = ^KG-^y (5.6) 

where \i m n is the number density of (m,n) strings, £„, „ is the energy of an individual (m,n) string, 
and g is the determinant of the metric. The source term 5? depends on the quantum numbers of the 
string gas, and the sum runs over all momentum and winding number vectors m and n, respectively 
(note that n and m are six-vectors, one component for each internal dimension). If the number of 
right-moving oscillator modes is given by N, then the source term for fixed m and n vectors is 

* = I(^) 2 " X>«*« + jhl + M+ 2 ( N ~ 0] • (5-7) 

To obtain this equation, we have made use of the mass spectrum of string states and of the level 
matching conditions. In the case of the bosonic superstring, the mass M of a string state with 
fixed m,n,N and N, where N and N are the number of right- and left-moving oscillator states, on a 
six-dimensional torus whose radii are given by b a is 

M 2 = { ^_^n 2 a b 2 a + 2(N + N-2), (5.8) 
°a a 

and the level matching condition reads 

N=(n,m)+N, (5.9) 

where (n,m) indicates the scalar product of n and m in the trivial metric. 

There are modes which are massless at the self-dual radius b a = 1. One such mode is the 
graviton with n = m = and N =\. The modes of interest to us are modes which contain winding 
and momentum, namely 

• N = \, (m,m) = 1, (m,n) = — 1 and (n,n) = 1; 

• N = 0, (m,m) = 1, (m,n) = 1 and (n,n) = 1; 

• N = (m,m) = 2, (m,n) = and (n,n) = 2. 
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Note that these modes survive in the Heterotic string theory, but do not survive the GSO [142] 
truncation in Type II string theories. 

In string theories which admit massless states (i.e. states which are massless at the self-dual 
radius), these states will dominate the initial partition function. The background dynamics will then 
also be dominated by these states. To understand the effect of these strings, consider the equation 
of motion ( |5.6| ) with the source term (fT7|). The first two terms in the source term correspond to an 
effective potential with a stable minimum at the self-dual radius. However, if the third term in the 
source does not vanish at the self-dual radius, it will lead to a positive potential which causes the 
radion to increase. Thus, a condition for the stabilization of b a at the self-dual radius is that the 



third term in ( p.7| ) vanishes at the self-dual radius. This is the case if and only if the string state is 
a massless mode. 

The massless modes have other nice features which are explored in detail in [ |160| ]. They 
act as radiation from the point of view of our three large dimensions and hence do not lead to a 
over-abundance problem. As our three spatial dimensions grow, the potential which confines the 
radion becomes shallower. However, rather surprisingly, it turns out the the potential remains steep 
enough to avoid fifth force constraints. 

In the presence of massless string states, the shape moduli also can be stabilized, at least in the 
simple toroidal backgrounds considered so far [ |164| ]. To study this issue, we consider a metric of 
the form 

ds 1 = dt 2 - dx 2 - G mn dy m dy n , (5.10) 

where the metric of the internal space (here for simplicity considered to be a two-dimensional 
torus) contains a shape modulus, the angle between the two cycles of the torus: 

G n =G 22 = \ (5.11) 

and 

G n = G 2l = sin0 , (5.12) 

where 6 = corresponds to a rectangular torus. The ratio between the two toroidal radii is a second 
shape modulus. However, from the discussion of the previous subsection we already know that each 
radion individually is stabilized at the self-dual radius. Thus, the shape modulus corresponding to 
the ratio of the toroidal radii is fixed, and the angle is the only shape modulus which is yet to be 
considered. 

Combining the 00 and the 12 Einstein equations, we obtain a harmonic oscillator equation for 
6 with 6 = as the stable fixed point. 

e + %K- l/2 e~ z *e = 0, (5.13) 

where A" is a constant whose value depends on the quantum numbers of the string gas. In the case 
of an expanding three-dimensional space we would have obtained an additional damping term in 
the above equation of motion. We thus conclude that the shape modulus is dynamically stabilized 
at a value which maximizes the area to circumference ratio. 

The only modulus which is not stabilized by string winding and momentum modes is the 



dilaton. Recently, it has been show [165] that a gaugino condensation mechanism (similar to those 
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used in string inflation model building) can be introduced which generates a stabilizing potential 
for the dilaton without interfering with the radion stabilization force provided by the string winding 
and momentum modes. 

In the next subsection we turn to the predictions of string gas cosmology for observations. 
These predictions do not depend on the details of the theory, but only on three inputs. The first 
is the existence of a quasi-static initial phase in thermal equilibrium. The second condition is the 
applicability of the Einstein field equations for fluctuations on infrared scales (scales of the order of 
1mm), many orders of magnitude larger than the string scale, and the third is a holographic scaling 
of the specific heat capacity. The role of the last condition will become manifest below. 



5.2 Fluctuations in String Gas Cosmology 
5.2.1 Overview 

At the outset of this section, let us recall the mechanism by which inflationary cosmology leads 
to the possibility of a causal generation mechanism for cosmological fluctuations which yields an 
almost scale-invariant spectrum of perturbations. The space-time diagram of inflationary cosmol- 
ogy is sketched in Figure |[ 

During the period of inflation, the Hubble radius 

l H {t) = a - (5.14) 

is approximately constant. In contrast, the physical length of a fixed co-moving scale (labelled by k 
in the figure) is expanding exponentially. In this way, in inflationary cosmology scales which have 
microscopic sub-Hubble wavelengths at the beginning of inflation are red-shifted to become super- 
Hubble-scale fluctuations at the end of the period of inflation. After inflation, the Hubble radius 
increases linearly in time, faster than the physical wavelength corresponding to a fixed co-moving 
scale. Thus, scales re-enter the Hubble radius at late times. 

Since inflation red-shifts any classical fluctuations which might have been present at the be- 
ginning of the inflationary phase, fluctuations in inflationary cosmology are generated by quantum 
vacuum perturbations. The fluctuations begin in their quantum vacuum state at the onset of in- 
flation. Once the wavelength exceeds the Hubble radius, squeezing of the wave-function of the 



fluctuations sets in (see [34, 33]). This squeezing plus the de-coherence of the fluctuations due to 



the interaction between short and long wavelength modes generated by the intrinsic non-linearities 



in both the gravitational and matter sectors of the theory (see [ J71| , [70| ] for recent discussions of this 
aspect and references to previous work) lead to the classicalization of the fluctuations on super- 
Hubble scales. 

Let us now turn to the cosmological background of string gas cosmology represented in Figure 



13[ This string gas cosmology background yields the space-time diagram sketched in Figure [15 
As in Figure ||, the vertical axis is time and the horizontal axis denotes the physical distance. For 
times t < tR, we are in the static Hagedorn phase and the Hubble radius is infinite. For t > the 
Einstein frame Hubble radius is expanding as in standard cosmology. The time is when the string 
winding modes begin to decay into string loops, and the scale factor starts to increase, leading to 
the transition to the radiation phase of standard cosmology. 
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Figure 15: Space-time diagram (sketch) showing the evolution of fixed co-moving scales in string gas 
cosmology. The vertical axis is time, the horizontal axis is physical distance. The solid curve represents the 
Einstein frame Hubble radius H 1 which shrinks abruptly to a micro-physical scale at tR and then increases 
linearly in time for t > tR. Fixed co-moving scales (the dotted lines labeled by k\ and h£) which are currently 
probed in cosmological observations have wavelengths which are smaller than the Hubble radius before tR. 
They exit the Hubble radius at times tj(k) just prior to tR, and propagate with a wavelength larger than the 
Hubble radius until they reenter the Hubble radius at times tf(k). 



Let us now compare the evolution of the physical wavelength corresponding to a fixed co- 
moving scale with that of the Einstein frame Hubble radius H~ l {t). The evolution of scales in 
string gas cosmology is identical to the evolution in standard and in inflationary cosmology for 
t > tR. If we follow the physical wavelength of the co-moving scale which corresponds to the 
current Hubble radius back to the time tR, then - taking the Hagedorn temperature to be of the order 
10 16 GeV - we obtain a length of about 1 mm. Compared to the string scale and the Planck scale, 
this is in the far infrared. 

The physical wavelength is constant in the Hagedorn phase since space is static. But, as we 
enter the Hagedorn phase going back in time, the Hubble radius diverges to infinity. Hence, as 
in the case of inflationary cosmology, fluctuation modes begin sub-Hubble during the Hagedorn 
phase, and thus a causal generation mechanism for fluctuations is possible. 

However, the physics of the generation mechanism is very different. In the case of inflationary 
cosmology, fluctuations are assumed to start as quantum vacuum perturbations because classical 
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inhomogeneities are red-shifting. In contrast, in the Hagedorn phase of string gas cosmology there 
is no red-shifting of classical matter. Hence, it is the fluctuations in the classical matter which 
dominate. Since classical matter is a string gas, the dominant fluctuations are string thermodynamic 
fluctuations. 



Our proposal for string gas structure formation is the following [ pip (see [ |16q ] for a more 
detailed description). For a fixed co-moving scale with wavenumber k we compute the matter 
fluctuations while the scale in sub-Hubble (and therefore gravitational effects are sub-dominant). 
When the scale exits the Hubble radius at time t\ (k) we use the gravitational constraint equations to 
determine the induced metric fluctuations, which are then propagated to late times using the usual 
equations of gravitational perturbation theory. Since the scales we are interested in are in the far 
infrared, we use the Einstein constraint equations for fluctuations. 

5.2.2 Spectrum of Cosmological Fluctuations 

We write the metric including cosmological perturbations (scalar metric fluctuations) and grav- 
itational waves in the following form (using conformal time 17) 

ds 2 = a 2 (ri) {(I +2<t>)dri 2 - [(I - 2®)8ij + h ij ]dx i dx j } . (5.15) 

We have fixed the gauge (i.e. coordinate) freedom for the scalar metric perturbations by adopting 
the longitudinal gauge and we have taken matter to be free of anisotropic stress. The spatial tensor 
hij(x,t) is transverse and traceless and represents the gravitational waves. 

Note that in contrast to the case of slow-roll inflation, scalar metric fluctuations and gravita- 
tional waves are generated by matter at the same order in cosmological perturbation theory. This 
could lead to the expectation that the amplitude of gravitational waves in string gas cosmology 
should be genetically larger than in inflationary cosmology. This expectation, however, is not real- 



ized [ |167Q since there is a different mechanism which suppresses the gravitational waves relative to 
the density perturbations (namely the fact that the gravitational wave amplitude is set by the ampli- 
tude of the pressure, and the pressure is suppressed relative to the energy density in the Hagedorn 
phase). 

Assuming that the fluctuations are described by the perturbed Einstein equations (they are 



not if the dilaton is not fixed [ 168| , |169| ]), then the spectra of cosmological perturbations <I> and 



gravitational waves h are given by the energy-momentum fluctuations in the following way [ ]166| ] 

(|<D(fc)| 2 ) = l6n 2 G 2 k- 4 (ST° (k)ST° (k)) , (5.16) 
where the pointed brackets indicate expectation values, and 

(\h(k)\ 2 ) = \6K 2 G 2 k-\8r j(k)Sr j(k)) , (5.17) 



where on the right hand side of ( p.l7| ) we mean the average over the correlation functions with 
i 7^ j, and h is the amplitude of the gravitational waves 14 . 

Let us now use ( 5.16D to determine the spectrum of scalar metric fluctuations. We first calculate 
the root mean square energy density fluctuations in a sphere of radius R = k~ l . For a system in 



I4 The gravitational wave tensor hn can be written as the amplitude h multiplied by a constant polarization tensor. 
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thermal equilibrium they are given by the specific heat capacity Cy via 

(8p 2 ) = ^C v . (5.18) 
The specific heat of a gas of closed strings on a torus of radius R can be derived from the partition 



function of a gas of closed strings. This computation was carried out in [17C] with the result 



R 2 I 3 

C v 2——t—— . (5.19) 
T(\-T/T H ) 

The specific heat capacity scales holographically with the size of the box. This result follows 
rigorously from evaluating the suing partition function in the Hagedorn phase. The result, however, 
can also be understood heuristically: in the Hagedorn phase the string winding modes are crucial. 
These modes look like point particles in one less spatial dimension. Hence, we expect the specific 
heat capacity to scale like in the case of point particles in one less dimension of space 15 . 

With these results, the power spectrum P(k) of scalar metric fluctuations can be evaluated as 
follows 

P*(k) = ^k 3 \<i>(k)\ 2 (5.20) 

= SG 2 k- 1 < \8p(k)\ 2 > . 
= SG 2 k 2 < (8M) 2 > R 
= 8G 2 k~ 4 < (dp) 2 > R 



8G 



T 1 

1 1 - T/T H 



where in the first step we have used ( 5. 16[ ) to replace the expectation value of | <!>(&) | 2 in terms of 



the correlation function of the energy density, and in the second step we have made the transition 
to position space 

The first conclusion from the result ( p. 20 ) is that the spectrum is approximately scale-invariant 



(P(k) is independent of k). It is the 'holographic' scaling Cy(R) ~ R 2 which is responsible for the 
overall scale-invariance of the spectrum of cosmological perturbations. However, there is a small k 
dependence which comes from the fact that in the above equation for a scale k the temperature T 
is to be evaluated at the time tj(k). Thus, the factor (1 — T /Tu) in the denominator is responsible 
for giving the spectrum a slight dependence on k. Since the temperature slightly decreases as time 
increases at the end of the Hagedorn phase, shorter wavelengths for which tj(k) occurs later obtain 
a smaller amplitude. Thus, the spectrum has a slight red tilt. 

5.2.3 Key Prediction of String Gas Cosmology 



As discovered in [167], the spectrum of gravitational waves is also nearly scale invariant. 
However, in the expression for the spectrum of gravitational waves the factor (\ — T/T H ) appears 
in the numerator, thus leading to a slight blue tilt in the spectrum. This is a prediction with which 
the cosmological effects of string gas cosmology can be distinguished from those of inflationary 



I5 We emphasize that it was important for us to have compact spatial dimensions in order to obtain the winding modes 
which are crucial to get the holographic scaling of the thermodynamic quantities. 
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cosmology, where quite generically a slight red tilt for gravitational waves results. The physical 
reason is that large scales exit the Hubble radius earlier when the pressure and hence also the 
off-diagonal spatial components of T^ v are closer to zero. 

Let us analyze this issue in a bit more detail and estimate the dimensionless power spectrum 
of gravitational waves. First, we make some general comments. In slow-roll inflation, to leading 
order in perturbation theory matter fluctuations do not couple to tensor modes. This is due to the 
fact that the matter background field is slowly evolving in time and the leading order gravitational 
fluctuations are linear in the matter fluctuations. In our case, the background is not evolving (at 
least at the level of our computations), and hence the dominant metric fluctuations are quadratic in 
the matter field fluctuations. At this level, matter fluctuations induce both scalar and tensor metric 
fluctuations. Based on this consideration we might expect that in our string gas cosmology scenario, 
the ratio of tensor to scalar metric fluctuations will be larger than in simple slow-roll inflationary 
models. However, since the amplitude h of the gravitational waves is proportional to the pressure, 
and the pressure is suppressed in the Hagedorn phase, the amplitude of the gravitational waves will 
also be small in string gas cosmology. 

The method for calculating the spectrum of gravitational waves is similar to the procedure out- 
lined in the last section for scalar metric fluctuations. For a mode with fixed co-moving wavenum- 
ber k, we compute the correlation function of the off-diagonal spatial elements of the string gas 



energy-momentum tensor at the time tt(k) when the mode exits the Hubble radius and use ( 5.17 ) to 
infer the amplitude of the power spectrum of gravitational waves at that time. We then follow the 
evolution of the gravitational wave power spectrum on super-Hubble scales for t > tj(k) using the 
equations of general relativistic perturbation theory. 



The power spectrum of the tensor modes is given by ( p.l7| ): 



P h (k) = 167r 2 G 2 rt J (Sr' j(k)Sr } {k)) 



(5.21) 



for i ^ j. Note that the k 3 factor multiplying the momentum space correlation function of T'j gives 
the position space correlation function C'/j(R) , namely the root mean square fluctuation of T'j in 
a region of radius R = k~ l (the reader who is skeptical about this point is invited to check that the 
dimensions work out properly). Thus, 



P h (k) 



l6n 2 G 2 k- 4 C/j(R) 



(5-22) 



The correlation function C'fj on the right hand side of the above equation follows from the thermal 
closed string partition function and was computed explicitly in [171] (see also [|ll|]). We obtain 



P h (k) ~ \6x 2 G 2 -^{\-T/T H )ln 2 



[1-T/Th) 



(5.23) 



which, for temperatures close to the Hagedorn value reduces to 

4 



Ph(k) ~ [ l f] {l-T/T H )\n 2 



jSp(l-T/T H ) 



(5-24) 



This shows that the spectrum of tensor modes is - to a first approximation, namely neglecting the 
logarithmic factor and neglecting the ^-dependence of T(ti(k)) - scale-invariant. 
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On super-Hubble scales, the amplitude h of the gravitational waves is constant. The wave 
oscillations freeze out when the wavelength of the mode crosses the Hubble radius. As in the 
case of scalar metric fluctuations, the waves are squeezed. Whereas the wave amplitude remains 
constant, its time derivative decreases. Another way to see this squeezing is to change variables to 

V0?,x) = a(ri)h(ri,x) (5.25) 
in terms of which the action has canonical kinetic term. The action in terms of y becomes 

S = 5 / rf4 * (V 2 - W,iW,i + ^V 2 ) (5-26) 
from which it immediately follows that on super-Hubble scales y ~ a. This is the squeezing of 



gravitational waves Q72Q. 

Since there is no /^-dependence in the squeezing factor, the scale-invariance of the spectrum at 
the end of the Hagedorn phase will lead to a scale-invariance of the spectrum at late times. 

Note that in the case of string gas cosmology, the squeezing factor z{t\ ) does not differ substan- 
tially from the squeezing factor a(r\) for gravitational waves. In the case of inflationary cosmology, 
z(t]) and a(r\) differ greatly during reheating, leading to a much larger squeezing for scalar metric 
fluctuations, and hence to a suppressed tensor to scalar ratio of fluctuations. In the case of string 
gas cosmology, there is no difference in squeezing between the scalar and the tensor modes. 

Let us return to the discussion of the spectrum of gravitational waves. The result for the power 



spectrum is given in (5.24), and we mentioned that to a first approximation this corresponds to 



a scale-invariant spectrum. As realized in [167], the logarithmic term and the ^-dependence of 



T(ti(k)) both lead to a small blue-tilt of the spectrum. This feature is characteristic of our scenario 
and cannot be reproduced in inflationary models. In inflationary models, the amplitude of the 
gravitational waves is set by the Hubble constant H. The Hubble constant cannot increase during 
inflation, and hence no blue tilt of the gravitational wave spectrum is possible. 

A heuristic way of understanding the origin of the slight blue tilt in the spectrum of tensor 
modes is as follows. The closer we get to the Hagedorn temperature, the more the thermal bath is 
dominated by long string states, and thus the smaller the pressure will be compared to the pressure 
of a pure radiation bath. Since the pressure terms (strictly speaking the anisotropic pressure terms) 
in the energy-momentum tensor are responsible for the tensor modes, we conclude that the smaller 
the value of the wavenumber k (and thus the higher the temperature r (?,•(&)) when the mode exits 
the Hubble radius, the lower the amplitude of the tensor modes. In contrast, the scalar modes are 
determined by the energy density, which increases at r (?,(&)) as k decreases, leading to a slight red 
tilt. 

The reader may ask about the predictions of string gas cosmology for non-Gaussianities. The 



answer is [ |172[ ] that the non-Gaussianities from the thermal string gas perturbations are Poisson- 



suppressed on scales larger than the thermal wavelength in the Hagedorn phase. However, if the 



spatial sections are initially large, then it is possible that a network of cosmic superstrings [173] 



will be left behind. These strings - if stable - would achieve a scaling solution (constant number 



of strings crossing each Hubble volume at each time [86, 87, pa]). Such strings give rise to linear 



discontinuities in the CMB temperature maps [174], lines which can be searched for using edge 



detection algorithms such as the Canny algorithm (see [175] for recent feasibility studies) 
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5.3 Problems of String Gas Cosmology 

The main problem of the current implementation of string gas cosmology is that it does not 
provide a quantitative model for the Hagedorn phase. The criticisms of string gas cosmology raised 



in [1169J, 176] are based on assuming that dilaton gravity should be the background. However, as 



mentioned earlier, this is not a reasonable assumption since the coupling of string gas matter to 
dilaton gravity is not S -duality invariant. The quasi-static early Hagedorn phase should also have 
constant dilaton. In addition, we do not expect that at high densities such as the Hagedorn density 
an simple effective action such as that of dilaton gravity will apply. Nevertheless, to make the 
string gas cosmology scenario into a real theory, it is crucial to obtain a good understanding of the 
background dynamics. For some initial steps in this direction see [ 149]. 

A second problem of string gas cosmology is the size problem (and the related entropy prob- 
lem). If the string scale is about 10 17 GeV as is preferred in early heterotic superstring models, then 
the radius of the universe during the Hagedorn phase must be many orders of magnitude larger than 
the string scale. Without embedding string gas cosmology into a bouncing cosmology, it seems un- 
natural to demand such a large initial size. This problem disappears if the Hagedorn phase is 
preceded by a phase of contraction, as in the model of [177]. In this case, however, it is non-trivial 
to arrange that the Hagedorn phase lasts sufficiently long to maintain thermal equilibrium over the 
required range of scales. 

It should be noted, however, that some of the conceptual problems of inflationary cosmology 
such as the trans-Planckian problem for fluctuations, do not arise in string gas cosmology. As in 
the case of the matter bounce scenario, the basic mechanism of the scenario is insensitive to what 
sets the cosmological constant to its observed very small value. 



6. Conclusions 

In these lectures I have argued that it is possible to explore physics of the very early universe 
using current cosmological observations. Given the large amount of new data which is expected 
over the next decade, early universe cosmology will be a vibrant field of research. The information 
about the very early universe is transferred to the current time mainly via imprints on the spectrum 
of cosmological fluctuations. Thus, the theory of cosmological fluctuations is the key tool of mod- 
ern cosmology. It links current observations with early universe cosmology. The theory can be 
applied to any background cosmology, not just inflationary cosmology. I have illustrated the appli- 
cation of the theory of cosmological perturbations to three scenarios of the very early universe, the 
inflationary scenario, the matter bounce paradigm, and string gas cosmology. 

The Hubble radius plays a key role in the evolution of fluctuations. On sub-Hubble scales the 
microphysical forces dominate, whereas on super-Hubble scales matter forces freeze out and grav- 
ity dominates. In order to have a causal mechanism of structure formation, it is therefore crucial 
that scales of current interest in cosmology originate at very early times inside the Hubble radius, 
and that they then propagate over an extended period of time on super-Hubble scales. This propaga- 
tion on super-Hubble scales is required in order to obtain the squeezing of the perturbations which 
is required to explain the acoustic oscillations in the angular power spectrum of CMB anisotropies. 

The inflationary scenario is the current paradigm of early universe cosmology. It explains 
important conceptual problems of Standard Big Bang Cosmology, and most importantly it provided 
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a mechanism to explain the origin of structure in the universe via causal physics. The perturbations 
begin as quantum vacuum fluctuations on sub-Hubble scales. They exit the Hubble radius during 
the phase of inflationary expansion of space and are then squeezed. Inflationary cosmology has 
been predictive: it predicted the detailed shape of the angular power spectrum of CMB anisotropics 
more than 15 years before the observations. However, it is important to realize that any theory 
which provides a scale-invariant spectrum of primordial curvature fluctuations and which admits 
a period during which the scales which are observed today evolve outside of the Hubble radius 
will have the same predictions, as already realized more than a decade before the development 



of inflationary cosmology by Sunyaev and Zel'dovich [ ]178| ] and by Peebles and Yu [ |179| ] (see 
also [36, ^]). Thus, it is incorrect to claim that the current observations confirm the inflationary 
scenario. In fact, current realizations of inflationary cosmology suffer from some basic conceptual 
problems. This motivates the search for an improved understanding of the very early universe. 

I have described two alternative early universe scenarios, the "matter bounce", and "string gas 
cosmology". Both lead to a scale-invariant spectrum of curvature fluctuations and involve squeez- 
ing of the perturbations on super-Hubble scales and are thus in agreement with current observations. 
The challenge for alternative scenarios is to identify clean predictions with which the new scenar- 
ios can be distinguished from those of inflation. In the case of the matter bounce the prediction we 
have identified is a particular shape of the bispectrum, in the case of string gas cosmology the key 
prediction is a slight blue tilt in the spectrum of gravitational waves. 

It is important to realize that physics of the very early universe is determined by physics of 
the highest energy scales, scales many orders of magnitude higher than those which are probed 
in accelerator experiments. Thus, a better understanding of quantum gravity (e.g. string theory) 
might lead to completely new possibilities for cosmology. In turn, it is only through cosmology 
that physics of the highest energies can be probed. 
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